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ABSTRACT 


This report documents the continued investigation of 
the design of Ride Quality Augmentation Systems ( RQAS ) 
for commuter aircraft. The purpose of these RQAS is the 
reduction of the vertical and lateral acceleration 
response of the aircraft due to atmospheric turbulence 
by the application of active control. The current in- 
vestigations include the refinement of the sample data 
feedback control laws based on the control-rate-weight- 
ing and output-weighting optimal control design tech- 
niques. These control designs were evaluated using 
aircraft time simulations driven by Dryden spectra tur- 
bulence. Fixed gain controllers were tested throughout 
the aircraft operating envelope. The preliminary design 
of the hardware modifications necessary to implement and 
test the RQAS on a commuter aircraft is included. These 
include a separate surface elevator and the flap modifi- 
cations to provide both direct lift and roll control. 

The results indicate that vertical acceleration 
reductions of 45% and lateral reductions of more than 
50% are possible. A fixed gain controller appears to be 
feasible with only minor response degradation. 
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NOTATION 


For this report, bold-face upper-case letters are used 
to denote matrices while bold-face lower-case letters 
are used to denote vectors. The prime symbol (') is used 
with matrices to denote matrix transpose and the 
superscript (-1) is used to denote matrix inverse. A dot 
over a variable is used to denote differentiation with 
respect to time. 
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CHAPTER 1. 


INTRODUCTION 


Since the 1978 federal deregulation of the major 
air carriers, there has been an expansion in the small- 
er, commuter class air carriers into the routes that are 
not profitable for the larger carriers. With the re- 
newed market for small (15-50 passengers) aircraft, 
comes renewed interest in technological advances for 
small aircraft. While many new advances are being in- 
corporated into existing aircraft and new designs, one 
area has received little attention, that of ride smooth- 
ness or ride quality. Due to the inherent characteris- 
tics of smaller aircraft, they are more susceptible to 
atmospheric gusts. This report gives an update on the 
work on a Ride Quality Augmentation System (RQAS) by the 
University of Kansas Flight Research Laboratory (KU- 
FRL) . RQAS is the implementation of an active digital 
flight control system for the expressed purpose of re- 
ducing aircraft vertical and lateral accelerations due 
to atmospheric turbulence. All of the research on the 
RQAS conducted at KU-FRL was done under the support and 
guidance of NASA Langley Research Center. 

The initial investigations which led to the current 
work involved a study of previous ride quality research 
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and a feasibility study to determine the best approach 
to implementing a active digital control system. (Refer- 
ence 1) The first phase of the current work (Reference 
2) began the theoretical design phase of a digital 
controller to be implemented on a Cessna 402B aircraft. 
This phase included the development of the Interactive 
Control Augmentation Design (ICAD) program which incor- 
porates classical and optimal control design techniques 
along with several different analysis techniques into 
one package. Using the ICAD program, longitudinal RQAS 
controllers were designed and evaluated in batch simula- 
tions, on the KU-FRL hybrid simulator and on the NASA 
Langley Research Center nonlinear moving-base simulator. 

The current work, described in this report, in- 
volves the continuation of the work of phase 1. The 
lateral and longitudinal RQAS controllers are designed 
using both the control-rate-weighting and output-weight- 
ing sampled-data optimal control techniques. A basic 
experimental system is defined. Preliminary design of 
the necessary aircraft modifications for installing 
direct lift flaps and a separate surface elevator are 
described. These surfaces will provide the necessary 
active control power while preventing annoying feedback 
to the pilot. It was originally intended to include a 
separate surface rudder to provide control in the later- 
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al direction. All theoretical designs are based upon a 
separate surface rudder which is 33% as effective as the 
full rudder. This is the separate surface rudder re- 
fered to throughout the current work. Since the theore- 
tical work was completed, the decision was made to use 
the entire rudder for lateral control. Since all models 
used are linear, all results for the separate surface 
rudder can be converted directly to full rudder results 
by a simple scale factor. 

The basic RQAS design problem is presented in Chap- 
ter 2. This includes models and chosen design para- 
meters. The sampled data optimal control designs are 
given in Chapter 3. The fixed gain controllers are 
included in Chapter 3. Chapter 4 outlines the work 
completed on the system modifications including the 

preliminary separate surface designs. 

I 
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CHAPTER 2. 


PROBLEM DEFINITION 

A ride quality augmentation system ( RQAS ) is an 
active control system which improves passenger and 
flight crew comfort. This type of system is generally 
designed to suppress an aircraft’s rigid body response 
to moderate to heavy continuous atmospheric turbulence. 
Although no standard criteria now exist for predicting 
comfort, several mathematical models of passenger re- 
sponse to aircraft motion have been developed and all 
agree that the dominant factors are the vertical and 
lateral accelerations (Reference 2). 

In this chapter, the RQAS problem is defined. The 
small-perturbation equations of motion for the Cessna 
402B aircraft are presented and the system design 
criteria are established. The atmospheric turbulence 
model to be used in the control law analysis is also 
discussed. 

2.1 CESSNA 40 2B AIRCRAFT EQUATIONS OF MOTION 

The RQAS control law designs presented in this 
report were designed using variations of the optimal 
linear quadratic regulator. The fundamental assumption 
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in applying these techniques to RQAS design is that the 
aircraft dynamics can be described by a set of linear, 
small-perturbation equations of motion in a state-matrix 
form; 

x = Ax + Bu , (2.1) 

y » Cx + Du (2.2) 

where x is the aircraft state vector, u is the aircraft 
control vector and y is the aircraft output vector. The 
matrices A, B, C and D are constant coefficient 
matrices . 

The linear, small-perturbation mathematical models 
used in this study were furnished for the Cessna 402B 
aircraft (Figure 2.1) by NASA Langley Research Center 
(LaRC) . These models were obtained from a nonlinear 
simulation model using standard NASA LaRC techniques 
(Reference 3). The primary assumptions which restrict 
the validity of these models are: 

1. The airframe is a rigid body: 

2. The earth is an inertial reference frame: 

3. The aircraft mass and mass distribution are 
constant : 

4. The XZ-plane is a plane of symmetry: 

5. The flow is quasi-steady: 

6. The effect of engine gyroscopics is negligible: 
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CLEARANCE IS 8.15 
INCHES. 



2. TOTAL WING AREA. 
INCLUDING NACELLES 
ANO FUSELAGE WITHIN 
THE WING PLANFORM. 


IS 195.72 SQUARE FEET. 


3. MINIMUM TURNING 


RADIUS IS 54'Q”. 



Three view 




7. The steady-state conditions are for straight 
line trimmed flight: 

8. The perturbations from steady-state are small: 

9. The equations are written with respect to the 
stability axis system. 

The state-matrix models were furnished for a 
coupled six-degree-of-f reedom linear model. These equa- 
tions were decoupled into the standard longitudinal and 
lateral-directional modes by assuming that the long- 
itudinal forces and moments due to lateral perturbations 
are negligible and vice versa. This assumption is valid 
since the coupling terms in the provided models were 
indeed small and the eigenvalues of the decoupled mat- 
rices matched those of the fully coupled matrices. 

The derivation of the equations of motion can be 
found in any standard text on airplane flight mechanics 
(Reference 4). The conversion of these equations to 
state matrix form is discussed in Appendix A. 

The standard controls for the 402B are the elevator 
for pitch control, the rudder for yaw control and the 
ailerons for roll control. The longitudinal controls 
used by the RQAS are direct lift flaps and a separate 
surface elevator. The lateral-directional controls used 
by the system are a separate surface rudder and differ- 
ential flaps. The outboard flaps are used as both longi- 
tudinal and lateral controls. Direct lift is achieved 
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when both the right and left flaps operate in the same 
direction and roll moments are developed by operating 
them differentially. The separate surface controls are 
provided so that the augmentation system is completly 
separate from the aircraft primary control system. A de- 
tailed description of how these controls are to be 
implemented on a Cessna 402B is contained in Chapter 4. 

The longitudinal equations of motion can therefore 
be written in the form of Equation (2.1) as follows: 


Za 

Zu 

Z q 

Z6 

Xi 


0 

x$ 

Mi 

Mi 

Mq 

Mi 


0 

1 

0 


a 


Z ‘se 

Z *f‘ 

u 

+ 

X *se 

n f 

q 


MJ 

se 

M4 f 

0 


0 

0 
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se 



(2.3) 


where the elements of the longitudinal A and B matrices 
are related to the standard aircraft longitudinal 
dimensional derivatives (Appendix A). 

The lateral-directional equations of motion can 
also be written in the form of Equation (2.1) as fol- 
lows: 
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Y & 

y p 

Y f 

n 

L b 

l p 

Lf 

0 

H 

Np 

Mf 

0 

0 

1 

tanQi 

0 


6 


Y *af 

Y *sr 

P 

+ 

L ‘dt 

L *sr 

r 


N *df 

N *sr 

♦ 


0 

0 


*df 

«sr 


(2.4) 


where the elements of the lateral-directional A and B 
matrices are related to the standard aircraft lateral- 
directional dimensional stability derivatives (Appendix 
A) . 


2.2 OUTPUT EQUATION FORMULATION 

As stated earlier, the dominant factors which in- 
fluence aircraft passenger comfort are the vertical and 
lateral accelerations. The three angular rates, pitch 
rate q, roll rate p, and yaw rate r, are also factors to 
a smaller extent. The yaw rate and roll rate can be 
especially uncomfortable in combination. Therefore, 
while the primary objective of a ride quality augmenta- 
tion system is to reduce the gust induced vertical and 
lateral accelerations, the systems designed should also 
try to maintain at least the open loop state response. 
With this in mind, the output vector is simply the state 
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vector with the addition of the two accelerations. The 
designer therefore has direct control over all of the 
pertinent (from a ride quality point of view) motion 
variables. 

The aircraft longitudinal state, control, and 
output vectors are: 

x' = (a, u, q, 8) 

u ' = ( 6 S e ’ ® f ) 

y' = (a z , a, u, q, 6) 

and the aircraft lateral-directional state, control, 
and output vectors are: 

x' = (B, p, r, <t>) 

u' = Udf* 5 sr ) 

y ' = (ay, B, p, r , <t> ) . 

The aircraft output equations must be written in the 
form of Equation (2.2). To do this, the accelerations a z 
and a y must be written as linear combinations of the 
states and controls. This is done by manipulating the 
following equations of motion from Appendix A: 


m(w - Uiq) = -mgysinyi + f a + ^T_ 

z z 

m(v + Uir) = mg4>cos0i + f^ + f<r 


(2.5) 

( 2 . 6 ) 


10 


can be 


The perturbed vertical acceleration a z 
obtained from Equation (2.5) by introducing the 
expressions; 

• • 

w = Uia and y * a - 0 : 

m(U]a - Uiq - gasiny! + gBsinyi) ® f A + f T 

z z 

Solving for a z results in: 

a z = ( f A, + ft )/m = u^a - gasinyi - Uiq + gesinyi 

Z (2.7) 

Introducing the expression for a from Equation (2.3) 

; into Equation (2.7) leads to the following equation for 

a z i 

a Z = Za'<* + ' u + Zq ' q + zy& + Z$'6 se + Z$’<5 f . (2.8) 

A similar derivation can be used to find the 
expression for ay. Substituting the equation v = U^B 
into Equation (2.6) and rearranging results in the 
equation 

i ay = (f& + f T )/m = u^B + Uir - g<|>cos0i (2.9) 

y y 

i Eliminating B from this expression using Equation (2.4) 

leads to the following equation for ay? 
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a y = Yfc'B + Yp ' p + Yyr + Y<J> ’ 4> + Y$^ d f + Y^5 sr . 


(2.10) 


The definitions of the modified dimensional stability 
derivatives contained in Equations (2.8) and (2.10) are 
summarized in Table 2.1. 


Table 2.1 Modified Dimensional Stability Derivatives 
For Use In The Aircraft Output Equations 

Lateral -Direct ional 


Longitudinal 

Z<$t ' = UiZi - gsinyi 

Zu' = UiZ{, 

Z^' = (Z^ - l)Ui 
zy = UiZ§ + gsinyi 

zy = u : z^ 

se se 

z$’ = UiZ4 f 

Note: The primed values (') 


yy = uiYfe 

Yp" = UlYp 

Yf' = (Yf + l)Ui 

yy = UiY6 - gcosOi 

Yk ' = UiY^ 

°se se 

Y$’ = UiY4 f 

are defined in Appendix A. 
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The longitudinal output equations can now be writ- 
ten in the form of Equation (2.2) as follows: 



Similarly, the lateral-directional output equations 
can be written as: 



13 


As an example of the models used in the RQAS designs, 
the longitudinal and lateral-directional equations of 
motion for the 402B during a sea level climb are 
presented in Tables 2.2 and 2.3. Models for all of the 
flight conditions investigated are included in Appendix 
B. 


2.3 DYNAMIC MODEL OF ATMOSPHERIC TURBULENCE 

Atmospheric turbulence may be described as 
individual patches in which the flow-field as seen by an 
airplane is a continuous random vector process composed 
of a steady mean value with turbulent fluctuations 
superposed. Each patch is assumed random, homogeneous 
(statistical properties are the same at each point in 
the field) and isotropic (statistical properties are 
independent of the axis orientation) . The turbulence 
intensity ,a, is used to distinguish one patch of turbu- 
lence from another. The frequency spectrum of the turbu- 
lence in each patch is related to this intensity. A more 
detailed explanation of the important concepts of turbu- 
lence as applied to the aircraft problem can be found in 
Reference 5. 

The spectral form for the random continuous turbu- 
lence model chosen for this project is the Dryden 
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Table 2.2 Cessna 402B Longitudinal Mathematical Model 
For Sea Level Climb 


x = Ax + Bu 


• 





- 


'a] 

a 


-1.3325 

-0.0014 

0.9189 

-0.0120 


• 

u 


12.7885 

-0.0228 

0.0000 

-32.0688 


u 

• 

q 


-6.4781 

0.0023 

-8.1525 

0.0406 

1 

1 

! 

, 

q 

• 

e 

— 


0.0000 

0.0000 

1.0000 

0.0000 i 


0 

i 


+ 


-0.0389 

0.0000 

-4.6678 

0.0000 


~l 

-0.2595 

-5.9354 

1.4135 

0.0000 

J 


$se 


y = Cx + Du 


V 


_ -283. 5580 

-0.2857 

-17.1010 

0.0780 


— 

a 


1.0000 

0.0000 

0.0000 

0.0000 


a 

u 

= 

0.0000 

1.0000 

0.0000 

0.0000 


u 

q 


0.0000 

0.0000 

1.0000 

0.0000 j 


q 

e! 

L J 

e 


0.0000 

0.0000 

0.0000 

j 

l.OOOOj 



-8.2010 

-54.7130 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

o.oooo: 

j 

0.0000 

0.0000_, 


^se 
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Table 2.3 Cessna 402B Lateral-Directional 

Mathematical Model For Sea Level Climb 


x = Ax + Bu 


p- 

8 


-0.1879 

0.0874 

-0.9971 

0.1505 


• 

P 


-3.7107 

-2.6275 

0.3918 

-0.0070 


• 

r 


3.7138 

-0.2901 

-0.3503 

-0.0065 


* 

l 

0.0000 

1.0000 

0.1700 

0.0000 



1 

8 1 

i 

P! 

I 

r 

4>l 


{ = 

Cx + Du 

~a~ 

' 

8 


P 

= 

r 

1 

| 

!_*_ 

i 

| 


-39.6290 

18.4390 

1.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 


0.0000 

0.0162 


-2.6247 

0.3362 


$df 

-0.0611 

-0.7013 


$sr 




- - 

0.0000 

0.0000 



0.6125 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


0.0193 

0.0000 

0.0000 

0.0000 

1 . 0000 _ 

3.4133 

0.0000 

0.0000 

0.0000 

0.0000 


5 df 

«sr 
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spectral form. The Dryden spectral form has the advan- 
tage over the Von Karman form in that it is easy to 
implement in the time domain. The basic approach is to 
set up a random gust disturbance file whose points have 
a root mean square value of one and a frequency spectrum 
whose behavior approximates that of the Dryden spectrum. 
The reader is referred to Reference 2 for a detailed 
discussion on how atmospheric turbulence is implemented 
in the simulation portion of the ICAD program. 

2.4 RQAS DESIGN CRITERIA 

As stated earlier in this chapter, the main object- 
ive of a ride quality augmentation system is to reduce 
the vertical and lateral accelerations without signi- 
ficantly degrading the other state variable responses 
from the aircraft's open loop response. Also, the three 
angular rates, q, p, and r contribute to passenger 
comfort to some degree. Therefore, any reduction in 
these rates is also desirable. 

The two factors which limit the performance of an 
RQAS are the rate and deflection limits for the control 
surfaces. For the 402B aircraft, the flaps are limited 
to deflections of 15 degrees up and 45 degrees down. 
These limits allow the RQAS to deflect +15 degrees about 
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any trimmed flap setting. The flap rate limit is set at 
120 degrees/second. The separate surface elevator de- 
flection and rate limits are set at +5 degrees and 50 
degrees/second respectively and the rudder deflection 
and rate limits are set at +5 degrees and 50 degrees/se- 
cond respectively. These limits are representative of 
available actuators (Reference 6) and will be shown to 
be adequate to provide good RQAS performance (Reference 
1 ) . 

The RQAS design criteria outlined above are listed 
in Table 2.4. 
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Table 2.4 


RQAS Design Criteria 


Longitudinal Mode 

Variable 
a 2 ( rms ) 

$ f 
<S S e 
l«fl 
I ^sel 

a,u,q,0 


Criterion 

< 3.54 (ft/sec 2 ) 

< 120 (deg/sec) 

< 50 (deg/sec) 

< 15 (deg) 

< 5 (deg) 

As close to the open loop 
values as possible. 


Lateral-Directional Mode 


Variable 
a y ( rms) 

5 df 

$sr 

I «df I 

l«rl 

P»r 

Bf ♦ 


Criterion 

< 50% of open loop 

< 120 (deg/sec) 

< 50 (deg/sec) 

< 15 (deg) 

< 5 (deg) 

Any reduction is desirable. 

As close to the open loop 
values as possible. 
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2.5 DESIGN PARAMETERS AND FLIGHT CONDITIONS 


Based upon the preliminary results from phase 1 
(Reference 2), it was decided to use a digital, linear 
quadratic regulator formulation for the control system 
designs. The designs wil be evaluated using a digital 
simulation in both the time and frequency domains. The 
sample time for the control law designs is 0.02 seconds, 
with a computational delay time assumed to be 0.002 
seconds. Also servos with bandwidths of 10 radians/ 
second are assumed for each control. The gust intensity 
is selected by assuming a probability of exceedence of 
0.001. These design parameters are summarized in Table 
2.5. 

Five flight conditions were chosen to represent a 
complete mission of the 402B. Table 2.6 summarizes these 
flight conditions. 
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Table 2.5 RQAS Design Parameters 


Digital, Optimal Regulator Controller 
Sample Time: T s = 0.02 seconds 

Computational Delay Time: T^ = 0.002 seconds 


Control Servos: 


10 
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Table 2.6 Flight Conditions 


Aircraft: Cessna 402B (2 crew / 6 passenger commuter) 


Condition 

Speed 
( ft/sec) 

0.001 P.O.E 

°z = °y 
( f t/sec £ ) 

Takeoff at Sea Level 

184 

9.5 

Climb at Sea Level 

211 

9.5 

Climb at 5,000 Feet 

227 

10.2 

Cruise at 20,000 Feet 

358 

7.4 

Approach at Sea Level 

160 

9.5 
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CHAPTER 3. 


OPTIMAL CONTROL DESIGNS 

In phase 1 (Reference 2), the RQAS control laws 
were designed using the sampled-data control-rate- 
weighting design technique. This method allowed the 
vertical acceleration to be written as a state variable. 
In the current phase of this work, a second control 
structure, output-weighting, was used to investigate 
alternative designs. (Reference 7) Both the control- 
rate weighting and the output-weighting design tech- 
niques are implemented in the ICAD program. This chap- 
ter will summarize the equations upon which both output- 
weighting and control-rate-weighting are based. These 
techniques will then be applied to the Cessna 402B model 
for all flight conditions. Since the optimal gains vary 
with flight conditions, the topic of gain scheduling 
will be investigated. 

3.1 DESIGN TECHNIQUES 

While both control-rate-weighting and output- 
weighting allow the control designer to weight outputs 
which are linear combinations of states and controls, 
each technique has distinct advantages and disadvan- 
tages. Control-rate-weighting allows the designer to 
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weight the control surface deflection rates directly, 
while output-weighting only allows weights on the con- 
trol deflection and gives no direct access to the con- 
trol rates. Control-rate-weighting, however, has the 
disadvantage in that the surface positions must be 
measured and fed back. By contrast, only the states must 
be fed back in output-weighting. 

Using the linear, small-perturbation equations of 
motion in state-matrix form as given in Equations (2.1) 
and (2.2), the continuous and sampled-data standard 
optimal regulator problem can be defined. The equations 
necessary to solve the continuous and the sampled-data 
problem are summarized in Tables 3.1 and 3.2, respect- 
ively. Included in these tables are the state and 
output equations which defined the equations of motion 
for the aircraft, and the quadratic cost functional, 
which is the basis of the optimal control problem. This 
functional is minimized by defining the feedback gains, 
K, based upon the solution of the matrix Riccati equa- 
tion. 

To solve the control-rate-weighting optimal control 
problem, an extension of the standard optimal regulator 
problem is made. In addition to the matrices that 
weight the outputs and controls in the cost functional, 
Q and R, a third matrix is included in the functional. 
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The S matrix weights the control deflection rates. To 
use the ORACLS design routines (Reference 8) on this 
newly defined cost functional, some matrix manipulation 
is done. This manipulation creates several new matrices 
which are combinations of the previous matrices. This 
converts the problem into the standard optimal regulator 
problem for which a solution is known. The basic equa- 
tions and new matrix definitions are presented in Tables 
3.1 and 3.2 for the continuous and the sampled-data 
cases. For a complete derivation of the equations the 
reader is referred to Reference 2 or 7 . 

The output-weighting problem is solved by similar 
techniques. First, the matrices are manipulated to form 
the problem into one analogous to the standard regulator 
problem, and then the problem is solved by existing 
techniques. The matrices that result from this process 
are presented in Tables 3.1 and 3.2 also. Once again 
the reader is referred to Reference 7 for a complete 
discussion and derivation of the output-weighting prob- 
lem. 

Tables 3.1 and 3.2 obviously contain a vast amount 
of information, much more than has been discussed here. 
They are presented in such a condensed form for the 
purpose of completeness of the problem definition and 
solution. 
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In both formulations, the optimal controller de- 
signs require the solution of a nonlinear algebraic 
matrix equation. ICAD is structured to perform the 
matrix manipulations to construct the required matrices 
and then solve the nonlinear algebraic equation. 
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Table 3.2 Summary of Sampled Data Optimal Regulator 
Problems and Solutions (Reference 7) 








where 


<&(t,t n ) = exp[A(t-t n )] 


t 

r(t,t n ) = {j4>(t,x)dT}B 


T 

Q = /{*' (t,0)Q«(t,0)}dt 
0 

T 

R = RT + / { r ’ (t,0)Qr(t,0)}dt 
0 

T 

H = /{*' (t,0)Qr(t,0)}dt 
0 


Q S = C'QC 


T 

Qs - /{*' (t,0)Qs*(t f 0)}dt 

0 

T 

r = rt + J { r * (t,o)Q s r(t,o) }dt 
0 

T 

M = /{«' (t,0)Q s r(t,0)}dt 
0 

Table 3.2 Summary of Sampled Data Optimal Regulator 
Problems and Solutions (Continued) 
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and 


$ c (t) = exp(A c t) 

T 

r c (t) * J{* c (t) d t}B c 
0 


T 

Qc 3 

0 

T 

Re = R C T + J{ri(t)Q c r c (t)}dt 
0 

T 

M c = J{$i(t)Q c r c (t)}dt 
0 


T 

Qo = /{*’ (t)[C’QC]«(t)}dt 
0 

T t 

M q = /{$' (t) [C*QC]/$(T)B(T)dT + «’(t)C’QD}dt 
0 0 

T t t 

R 0 = (D ' QD+R) + J { JB ' (t)$' (T)dT(C'QC)/<KT)B(T)dT +H'+N}dt 

0 0 . 0 


t 

H * D'QCf$(T)B(T)dT 
0 


Table 3.2 Summary of Sampled Data Optimal Regulator 
Problems and Solutions (Concluded) 
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3.2 SUMMARY OF SAMPLED-DATA OPTIMAL CONTROL DESIGNS 

Using the techniques of Appendix A and Reference 3, 
the Cessna 402B equations of motion were modeled in 
state-space form for five flight conditions spanning the 
basic operating envelope of the aircraft. These condi- 
tions are 

1. Take-off at sea-level conditions 

2. Climb at sea-level 

3. Climb at 5000 feet 

4. Cruise at 20,000 feet 

5. Approach at sea-level. 

Both the longitudinal and lateral directional models for 
all five flight conditions are included in Appendix B. 
Using the ICAD program, the sampled-data control-rate- 
weighting and output-weighting design techniques pre- 
sented in Section 3.1 were applied to these linearized 
models. 

To illustrate the basic process and controller 
design, the optimal design of the RQAS for the sea-level 
climb flight condition will be presented in the next two 
sections. First, the control-rate-weighting designs 
will be presented, and then the output-weighting de- 
signs. A complete set of designs for all of the flight 
conditions is included in Appendix B. 
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3.2.1 CONTROL-RATE-WEIGHTING DESIGNS 


The RQAS performance using the sampled-data con- 
trol-rate-weighting design technique for the sea-level 
climb condition is presented. Performance is indicated 
by the maximum values and RMS values of the components 
of the output vector. Included are the open loop re- 
sponse, i.e. the response of the unaugmented aircraft to 
the gust field, and the response of the RQAS controlled 
aircraft. The system eigenvalues for both the unaugment- 
ed and augmented aircraft are given. The augmented 
aircraft eigenvalues are given in terms of W' -plane 
frequency and damping ratios. To evaluate the W' -plane 
eigenvalues, the controller is designed in the Z-plane. 
The Z-plane eigenvalues are then transformed into the 
W' -plane. The equivalence of the S-plane and the W'- 
plane is valid due to the small sample time used (Ts 
0.02 seconds). Also included are the final diagonal 
elements of the weighting matrices Q, R, and S, and the 
resulting optimal gain matrix K, used to implement the 
controller 

u = - K x 

where for the longitudinal motion 

x* = ( ot , u, q, 0, 6gg, 6f ) 

u ' = ( §se • ^ f ) 
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and for the lateral-directional motion 
x' = (8, Pr r, 4> , Sdf , S sr ) 

u ' = ( , 631- ) < 

3. 2. 1.1 LONGITUDINAL DESIGNS 

The results of the longitudinal control-rate- 
weighting designs, as presented in Table 3.3 for a sea- 
level climb, display many features common to the results 
for all of the longitudinal RQAS designs. First of all, 
there is a subtantial reduction in the vertical acceler- 
ations, both peak value and RMS value. That is to be 
expected since that is the expressed purpose of the 
RQAS. There is, however, a substantial increase in the 
other values when viewed as a percentage. When viewed in 
terms of absolute change in magnitude, the increase is 
not seen to be significant. For example, the maximum 
forward velocity perturbation, u, increases by 116% when 
the RQAS is engaged, but that increase is actually only 
a change of 3.2 feet per second, from 1.9 to 4.1 feet 
per second. The other values behave similarly. The 
time domain response portion of the table also indicates 
the control activity for the augmented aircraft. It can 
be seen that the maximum deflections for the flaps and 
separate surface elevator are at or below the standards 
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Open 

Loop 



Max. 

RMS 

Max. 

a z 

( f/s2 ) 

13.1 

6.3 

12.2 

a 

(deg) 

2.8 

1.3 

3.6 

u 

(f/s) 

1.9 

1.1 

4.1 

q 

(deg/s) 

1.8 

0.9 

2.7 

e 

(deg) 

1.6 

1.0 

2.0 

5 se 

(deg) 

- 

- 

4.3 

«f 

(deg) 

- 

- 

11.8 

5se 

(deg/s) 

- 

- 

11.4 

(deg/s) 

- 

- 

91.1 

* 
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nal RQAS Response 


ib — Control-Rate-Weighting 


seline 

IQAS 


Percent 

Reduction 


RMS 

Max. 

RMS 

3.4 

7% 

46% 

1.5 

-29% 

-15% 

2.4 

-116% 

-118% 

1.4 

-50% 

-56% 

1.1 

-25% 

-10% 

2.3 



__ 


3.4 


Table 3.3 Summary of Longitudinal RQAS Response 

for a Sea-Level Climb — Control-Rate-Weighting 
(Continued) 


System Eigenvalues : 


Open Loop 


Baseline RQAS 


Frequency Damping 


Frequency Damping 


Phugoid 

0.15 

0.04 

0.09 

0.79 

Phugoid 

0.15 

0.04 

0.09 

0.79 

Short Period 

2.4 

1.0 

6.0 

0.68 

Short Period 

7.1 

1.0 

6.0 

0.68 

Servo 

10.0 

1.0 

3.5 

1.0 

Servo 

10.0 

1.0 

11.8 

1.0 

Filter 

- 

- 

15.3 

1.0 

Filter 

- 

- 

63.8 

1.0 


Note 1: Baseline RQAS Frequency and Damping are based on 

roots of the W' -Plane. 


Note 2: Control-Rate-Weighting introduces 2 filter states. 


Matrix : 


Optimal 

K = 


Gain 
27.930 
171.39 


0.0123 

0.1719 


-2.5093 

14.655 


-13.634 

24.772 


8.8770 

3.0441 


6.1749 

38.537 


Weighting Matrices — Diagonal Elements : 


Q’ = (0.05, 10.0, 
R' = (0.12, 0.33) 
S’ = (0.02, 0.03) 


0.0001, 3.0, 20.0) 
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of Table 2.4. The maximum control rates are also well 
below the designated maximum rates. This is to be 
expected for in control-rate-weighting, as the name 
implies, the control rates can be directly weighted 
making maximum rate limits easy to enforce. Reference 7 
shows that the nature of the control-rate-weighting 
structure limits the maximum deflection rates. In that 
reference, the weighting of the rates was reduced to 
extremely small values but the actual deflection rates 
seem to approach an asymptotic value. This can be 
explained by the fact that the deflection rates are 
included in the cost functional, therefore, the de- 
flection cannot go to infinity while the cost func- 
tional, which penalizes control rates, is still mini- 
mized. 

The system eigenvalues are also presented in Table 
3.3 for the longitudinal mode. The eigenvalues are pre- 
sented in terms of frequencies and damping ratios. The 
open loop eigenvalues, i.e. the unaugmented aircraft 
eigenvalues, show that the aircraft phugoid mode is 
barely damped and the short period mode has separated 
into 2 real roots. The two servo eigenvalues are from 
the separate elevator and flap servos. The baseline 
RQAS augmented aircraft has very good characteristics. 
Both the phugoid and the short period have damping above 
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critical damping (0.707). There are two filter states 
introduced by control-rate-weighting in addition to the 
two servo states. 

3. 2. 1.2 LATERAL-DIRECTIONAL DESIGNS 

The lateral-directional RQAS design for the sea- 
level climb flight condition are representative of the 
characteristics of the lateral-directional designs for 
all five flight conditions. Good lateral-directional 
designs proved to be much easier to obtain than longi- 
tudinal designs. As can be seen in Table 3.4, there is 
over a 50% reduction in the RMS lateral acceleration 
when the baseline RQAS is implemented. Where this 
caused the remaining state values to increase in the 
longitudinal motion, all of the state values decrease in 
the lateral-directional motion. For example, the RMS 
value for the aircraft roll rate decreases by over 80%. 
The control surface deflections and deflection rates are 
also all well within the desired limits. These limits 
were actually the driving factor in the lateral-direc- 
tional motion. The lateral acceleration could be re- 
duced almost without limit if the deflections and rates 
were allowed to increase. The control-rate-weighting 
intrinsic rate limits, of course, held the control rates 
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in check. 

The lateral-directional system eigenvalues included 
in Table 3.4 indicate some of the reasons for the ease 
of design mentioned above. The unaugmented aircraft has 
a slightly unstable spiral mode and a lightly damped 
Dutch Roll mode. The instability of the spiral mode is 
not of much consequence as long as the eigenvalue is 
near the origin. The baseline RQAS augmented aircraft 
has considerably different eigenvalues which is expect- 
ed. The spiral has been substantially stabilized. The 
Dutch Roll mode has also been very well damped. There 
are also two filter states introduced. The excellent 
time domain response can be explained by the damping of 
the augmented Dutch Roll over the unaugmented Dutch 
Roll. 

A tabular summary of the baseline control-rate- 
weighting RQAS designs for all of the remaining flight 
conditions can be found in Appendix B. 


37 



Table 3.4 Summary of Lateral-Directional RQAS Response 

for a Sea-Level Climb — Control-Rate-Weighting 


Time Domain Response: 


Open 

Loop 


Baseline 

RQAS 


Percent 

Reduction 



Max. 

RMS 

Max. 

RMS 

Max. 

RMS 

a y 

(f/s2) 

6.1 

2.7 

3.5 

1.3 

43% 

52% 

B 

(deg) 

6.2 

2.8 

5.8 

1.9 

6% 

32% 

P 

(deg/s) 

6.4 

2.7 

1.5 

0.5 

77% 

81% 

r 

(deg/s) 

10.6 

4.8 

7.1 

2.8 

33% 

42% 

♦ 

(deg) 

3.2 

1.3 

1.0 

0.4 

69% 

69% 

«df 

(deg) 

- 

- 

4.8 

2.1 

— * 

" 

«sr 

(deg) 

- 

- 

14.7 

5.8 



«df 

(deg/s) 

- 

- 

30.0 

— 



5sr 

(deg/s) 

- 

- 

30.8 

- 

— 

— 
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Table 3.4 Summary of Lateral-Directional RQAS Response 

for a Sea-Level Climb — Control-Rate-Weighting 
(Continued) 


System Eigenvalues : 


Open Loop 


Baseline RQAS 


Frequency Damping Frequency Damping 


Spiral Mode 

0.02 

-1.0 

0.59 

1.0 

Roll Mode 

2.7 

1.0 

3.3 

1.0 

Dutch Roll 

2.1 

0.12 

2.8 

0.51 

Dutch Roll 

2.1 

0.12 

2.8 

0.51 

Servo 

10.0 

1.0 

5.5 

0.88 

Servo 

10.0 

1.0 

5.5 

0.88 

Filter 

- 

- 

12.6 

1.0 

Filter 

— 

- 

21.3 

1.0 

Note 1: Baseline RQAS 

roots of the 

Frequency 
W' -Plane. 

and Damping are 

based on 


Note 2: Control-Rate-Weighting introduces 2 filter states. 


Optimal Gain Matrix : 


K 


22.429 

-1.9138 


-9.4558 

3.7262 


-6.6832 

-14.873 


-12.139 

3.5634 


16.148 

-0.4145 


-0.6806 

7.2291 


Weighting Matrices — Diagonal Elements : 

Q' = (0.01, 15.0, 0.07, 0.5, 1.5) 

R’ = (2.007, 0.52) 

S’ = (0.007, 0.02) 
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3.2.2 OUTPUT-WEIGHTING DESIGNS 


The results of applying the sampled-data output- 
weighing design technique to the sea-level climb condi- 
tion are presented in Tables 3.5 and 3.6 for the lonitu- 
dinal and lateral-directional motions, respectively. 
These tables include all of the same items as the corre- 
sponding tables for the control-rate-weighting designs. 
The feedback gain matrix, K, along with the two weight- 
ing matrices, Q and R, are also included. There is no 
equivalent to the S weighting matrix in output-weight- 
ing. The state and control vectors for the longitudinal 
output- weighting designs are 
x' = (a, u, q, 0) 

U ' = ( 6 se r 5 f ) 

and for the lateral-directional design are 
x’ = (B, P, r, <j>) 
u' = ( Seif , S sr ) 

The state vector only has 4 elements for output-weight- 
ing as compared to 6 for control-rate-weighting. This 
difference can make output-weighting easier to implement 
because the control deflections are not part of the 
state vector and do not have to be measured. 
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3. 2. 2.1 LONGITUDINAL DESIGNS 


The results of the longitudinal output-weighting 
designs (Table 3.5) compare favorably to the control- 
rate-weighting designs. (Table 3.3) The RMS vertical 
acceleration is subtantially reduced over the unaugment— 
ed aircraft acceleration. The remaining values increase 
over the unaugmented aircraft values just as the con- 
trol-rate-weighting designs did. Again, the percent 
increases are rather large, but the magnitude to the 
results are not particular large. The control surface 
activity for the RQAS augmented aircraft also meets all 
of the design specifications. The activity of the out- 
put-weighting design is slightly greater than the acti- 
vity of the control-rate-weighting designs. The sepa- 
rate surface elevator is used more extensively, while 
the maximum control rates are both higher. This differ- 
ence can be attributed to two separate sources. First, 
the intrinsic rate control of control-rate-weighting 
holds those control rates down. Also, the design pro- 
cess is not unique, i.e. the designs both satisfy the 
design requirements, but the method by which the accel- 
eration reduction is achieved can vary considerably. 

The eigenvalues of the RQAS baseline output- 
weighting show good general characteristics. They are 
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Table 3.5 Summary of Longitudinal RQAS Response 

for a Sea-Level Climb — Output-Weighting 


Time Domain Response : 

Open 

Loop 


Baseline 

RQAS 


Percent 

Reduction 



Max. 

RMS 

Max. 

RMS 

Max. 

RMS 

a z 

(f/s2) 

13.1 

6.3 

11.8 

3.3 

10% 

48% 

a 

(deg) 

2.8 

1.3 

3.7 

1.6 

-32% 

-23% 

u 

( f/s) 

1.9 

1.1 

4.0 

2.3 

-111% 

-110% 

q 

(deg/s) 

1.8 

0.9 

2.5 

1.2 

-39% 

-33% 

6 

(deg) 

1.6 

1.0 

1.8 

1.1 

-13% 

-10% 

, 5 se 

(deg) 

- 

- 

5.0 

2.5 

— 

“ 

«f 

(deg) 

- 

- 

11.7 

5.4 

— 

“ 

6 se 

(deg/s) 

- 

- 

22.3 

— 

— 


h 

(deg/s) 

- 

- 

98.8 

— 

— 

— 
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Table 3.5 Summary of Longitudinal RQAS Response 

for a Sea-Level Climb — Output-Weighting 
(Continued) 


System Eigenvalues : 

Open Loop Baseline RQAS 

Frequency Damping Frequency Damping 


Phugoid 

0.15 

0.04 

0.09 

0.89 

Phugoid 

0.15 

0.04 

0.09 

0.89 

Short Period 

2.4 

1.0 

11.6 

0.66 

Short Period 

7.1 

1.0 

11.6 

0.66 

Servo 

10.0 

1.0 

1.9 

1.0 

Servo 

10.0 

1.0 

11.3 

1.0 


Note 1: Baseline RQAS Frequency and Damping are based on 

roots of the W' -Plane. 


Optimal Gain Matrix : 







K = 

-0.5014 

-0.0018 

-1.5609 

-2.5312 


4.5062 

0.0046 

0.5469 

0.9847 


Weighting Matrices — Diagonal Elements : 

Q’ = (0.05, 10.0, 0.0001, 18.0, 40.0) 
R' = (2.2, 1.0) 
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very stable and the oscillatory motion is well damped. 
The guaranteed stability is, of course, one of the 
benefits of optimal control. The phugoid mode is over 
critically damped and the short period mode is nearly 
so. The frequency of the short period is rather high. 
This is one of the characteristics observed in all of 
the output-weighting designs. The control-rate-weight- 
ing designs also have a high short period frequency, but 
not as high as the output-weighting designs. 


3. 2. 2. 2 LATERAL-DIRECTIONAL DESIGNS 


The results of the lateral -direct ional output- 
weighting designs of Table 3.6 have the same good char- 
acteristics as the control-rate-weighting designs. In 
addition to the marked decrease in the side acceler- 
ations, the other perturbation variables are also de- 
creased, sometimes by a large percentage. The RMS ac- 
celeration is reduced by 59%, which is much better than 
the 50% design criteria. This is achieved without any 
harm to the other variables. The only limit is the 
control surface deflection and rate limits. Since the 
control rates could not be directly weighted, the rate 
limits had to be enforced by using some engineering and 
increasing or decreasing the weights on some of the 
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Table 3.6 Summary of Lateral- 
for a Sea-Level Cli 




Max. 

RMS 

Max. 

a y 

(f/s2) 

6.1 

2.7 

3.4 

B 

(deg) 

6.2 

2.8 

6.0 

P 

(deg/s) 

6.4 

2.7 

1.8 

r 

(deg/s) 

10.6 

4.8 

6.8 

4> 

(deg) 

3.2 

1.3 

1.7 

6 df 

(deg) 

- 

- 

6.2 

«sr 

(deg) 

- 

- 

14.7 

, $df 
(deg/s) 

- 

- 

51.4 

®sr 

(deg/s) 

- 

- 

27.1 
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irectional RQAS Response 
b — Ouput-Weighting 


seline 

RQAS 


Percent 

Reduction 


RMS 

Max. 

RMS 

1.1 

44% 

59% 

1.9 

3% 

32% 

0.6 

72% 

78% 

2.7 

36% 

44% 

0.7 

47% 

46% 

2.8 

- 

- 

5.7 
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Table 3.6 Summary of Lateral-Directional RQAS Response 
for a Sea-Level Climb — Output-Weighting 
(Continued) 


System Eigenvalues: 


Open Loop 


Baseline RQAS 


Frequency Damping 


Frequency Damping 


Spiral Mode 

0.02 

-1.0 

0.78 

1.0 

Roll Mode 

2.7 

1.0 

7.4 

1.0 

Dutch Roll 

2.1 

0.12 

2.2 

0.58 

Dutch Roll 

2.1 

0.12 

2.2 

0.58 

Servo 

10.0 

1.0 

6.8 

0.87 

Servo 

10.0 

1.0 

6.8 

0.87 


Note 1: Baseline RQAS Frequency and Damping are based on 

roots of the W'-Plane. 


Optimal Gain Matrix : 



2.0004 

-O.i 

K = 




-0.0932 

1 . 


-0.6140 

-2.4060 


-1.0563 

1.1638 


Weighting Matrices — Diagonal Elements ; 

Q' * (0.05, 10.0, 0.1, 0.75, 10.0) 

R' = (7.0, 1.8) 
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other variables, most importantly the acceleration. The 
output-weighting designs did get slightly better accel- 
eration reductions, but the differential flap deflection 
rate was increased by about 60%. These factors must be 
balanced out in the final selection of an algorithm for 
implementation. 

The system eigenvalues of Table 3.6 are all well 
behaved. The Dutch Roll damping has been increased by a 
large margin over the unaugmented aircraft. The un- 
stable spiral mode has been strongly stabilized. A 
question might arise if this is too stable for the 
pilots. That is a subject recommended for future in- 
vestigations. The roll mode eigenvalue also has been 
pushed far into the left half of the W' -plane. This is 
comparable to the high frequencies developed in the 
short period mode by output-weighting. 

The output-weighting designs and performance sum- 
maries for all flight conditions can be found in Appen- 
dix B. 
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3.3 GAIN SCHEDULING 


The topic of gain scheduling for aircraft control 
systems is always a complex one. The need for gain 
scheduling arises from the variation of the aircraft 
open loop dynamics throughout its operation envelope. 
The feedback gains developed by optimal control tech- 
niques mirror the aircraft dynamics variation due to the 
dependence of the optimal regulator solution on the air- 
craft state equation. Therefore, many distinct sets of 
gains are needed to cover the flight envelope. There 
are two ways of solving this problem. One is to only 
use one set of gains for the entire flight envelope, and 
the other is to schedule the gains to coincide with the 
current flight condition. 

One set of fixed gains, constant for all flight 
conditions would be very easy to implement. The ques- 
tions that need to be addressed for fixed gain imple- 
mentations are related to the degradation of performance 
by the fixed gains over the optimal gains and the sta- 
bility of the fixed gain controlled aircraft in the 
entire flight envelope. If a fixed gain design is de- 
sired, the question of how to choose the gains arises. 
They may be chosen by averaging the gains for all of the 
available flight conditions, or, more likely, they may 
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be chosen by some sort of weighted average because all 
available designs may not be deemed equally important. 

The other method of handling the multiple sets of 
gains is gain scheduling. If gain scheduling is chosen, 
the scheduled gains must still be determined. Two major 
questions arise. What parameter (or parameters) should 
the gains be scheduled on and how should the gains be 
fit between design points? Common parameters used for 
gain scheduling include aircraft velocity, center of 
gravity location, and trimmed control surface position. 
The gains are usually fit between design points by some 
sort of curve fitting technique. This can introduce 
non-optimal gains though, i.e. the scheduled gains at 
the design points are not equal to the optimal gains. 
The questions of performance degradation and stability 
become important for these non-optimal gains. 

The baseline optimal designs for the Cessna 402B 
(Appendix B) were all optimized for a particular flight 
condition. Due to the variation of the 402B dynamics, 
the controller gains do vary "significantly" with flight 
condition. On the other hand, they are designs for a 
particular aircraft so the corresponding row and column 
positions in the gain matrices take on values that are 
in general of the same sign and order of magnitude 
although this is not universally true. All of these 
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items may be seen by comparing the optimal design gains 
in Tables B.21 through B.40 which cover all flight 
conditions for both the longitudinal and lateral-di- 
rectional controllers using control-rate-weighting and 
output-weighting . 

Initially it was thought that gain scheduling for 
the 402B would be necessary. The implementation prob- 
lems would be somewhat alleviated by the large amount of 
computing power available in the proposed onboard flight 
computer. If scheduled, the gains would be based upon a 
combination of forward velocity, trimmed flap de- 
flection, and center of gravity location. These para- 
meters would distinctly identify each basic flight con- 
dition allowing for appropriate gain scheduling. 

For the initial work on the RQAS gain scheduling 
problem, four sets of fixed gains were chosen and test- 
ed for feasibility. The four sets are 

1. Longitudinal Control-Rate-Weighting Fixed Gains 

2. Lateral-Direction Control-Rate-Weighting 
Fixed Gains 

3. Longitudinal Output-Weighting Fixed Gains 

4. Lateral-Directional Output-Weighting 
Fixed Gains. 

These fixed gains were determined by comparing all of 
the optimal gains for the five flight conditions and 
choosing a gain for each matrix position which would 
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come closest to most of the values. An average was not 
directly taken because the five flight conditions were 
not felt to be equally important. The take-off and 
approach conditions generally had optimal gains which 
were considerably different than the other three condi- 
tions. This was due to the aircraft forward velocity in 
these conditions and especially the trim flap deflection 
in the approach condition. 

The results of applying the fixed gain designs to 
the five flight conditions and the three longitudinal 
center of gravity locations follows in the next five 
sections. Included in these sections are figures which 
depict the accelerations and control surface activity 
for the baseline optimal designs and for the fixed gain 
designs. The fifth section investigates the stability 
of the fixed gain designs as compared to the baseline 
designs. 

The general results of this analysis indicate that 
no gain scheduling is needed. It is recommended that 
the issue be reexamined as refinements are made to the 
vehicle model. This would include models of the test 
aircraft with the modified controls. 
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3.3.1 LONGITUDINAL CONTROL-RATE-WE I GHT I NG FIXED GAINS 


The results of the longitudinal control-rate- 
weighting fixed gain tests are presented in Figures 3.1, 
3.2, and 3.3. Each plot contains the results of time 
simulation for three center of gravity locations. These 
are the quarter chord c.g., which is representative of 
the mid portion of the c.g. travel range, the forward 
most c.g. location, and the aft most c.g. location. All 
basic designs were done at the mid c.g. location. For 
the longitudinal designs, plots of the vertical RMS 
acceleration and the control surface activity, i.e. 
maximum deflections and rates, were made. 

The fixed gain matrix for the control-rate-weight- 
ing longitudinal design is 


28.0 

o 

• 

o 

-3.0 

-13.0 

9.0 

6.0 

170.0 

0.0 

14.5 

24.5 

3.0 

38.5 


Several points should be noticed. First, the column of 
zero elements corresponds to the velocity feedback col- 
umn. The actual elements existed in the optimal gains, 
but they were several orders of magnitude smaller than 
all of the other elements of the matrix. Therefore, 
they were set to zero. It might be argued that the 
velocity is several orders of magnitude larger than the 
other variables, but this velocity is a perturbation 
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around the steady state, so it should never be large, 
especially if the linear models are to hold. It is also 
obvious that all of the elements were rounded to "nice" 
values. This is true, but the method of determining 
these gains didn’t warrant any more accuracy. These 
gains could be directly compared to the baseline optimal 
gains. For some flight conditions, the gains are nearly 
identical, and for other conditions there are large 
differences. It remains to be seen how the differences 
affect the outcome. 

The RMS vertical acceleration is plotted in Figure 

3.1 for the fixed gain design, the baseline design and 
the unaugmented, open loop, aircraft. All results are 
generated with the same gust field so all results are 
directly comparable. The fixed gain designs perform 
nearly the same as the baseline designs. In the take- 
off flight condition, the fixed gain controller works 
better. The O.llg line is the RQAS upper design limit 
for acceptable vertical accelerations. The 0.055g line 
indicates a lower limit below which further reductions 
are deemed unnecessary and cost prohibitive. 

The excellent acceleration reduction of the fixed 
gain designs is not without its downfalls. These may be 
seen in the control surface activity shown in Figures 

3.2 and 3.3. All flight conditions except the take-off 
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□pen Loop 



(2-S/^H) 


Figure 3.1 Longitudinal Control-Rate-Weighting 
Fixed Gain Acceleration Performance 
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Flight Condition 



perform reasonably well. All values are near enough to 
the baseline to call them equal within the given accur- 
acy of the models. Maximum rates and deflections for 
both the flaps and the separate surface elevator in- 
crease for the take-off flight condition. That is the 
reason for the larger acceleration reduction. The 
rates, although larger, are still within the design 
limits. The surface deflections, though, are not within 
their limits. If the deflection limits were enforced, 
the vertical acceleration performance would deteriorate, 
but it is felt that it would not be too substantial. 
The same sort of argument could be made for the separate 
surface elevator deflection for the approach flight 
condition. 

The effect of center of gravity location is not 
particularly strong for the control-rate-weighting de- 
signs. The RMS accelerations are nearly identical for 
all of the flight conditions. The center of gravity 
location is a little more apparent on the surface acti- 
vity. The aft center of gravity location requires more 
flap deflection and rate than the other locations. The 
forward center of gravity location, on the other hand, 
uses the separate surface elevator more. 
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Maximum Flap Deflection (deg) 



Flight Condition 



12 3 4 

Flight Condition 


Figure 3.2 Longitudinal Control-Rate-Weighting 
Fixed Gain Direct Lift Flap Activity 
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Maximum Elevator Deflection (deg) Maximum Elevator Rate (deg/sec) 



Flight Condition 



1 2 3 4 5 

Flight Condition 


Figure 3.3 Longitudinal Control-Rate-Weighting Fixed 
Gain Separate Surface Elevator Activity 
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3.3.2 LATERAL-D I RECT I ONAL CONTROL-RATE- WE I GHT I NG 
FIXED GAINS 


The lateral-directional fixed gain design response 
can be seen in Figures 3.4 through 3.6. The control- 
rate-weighting fixed gains are 


K = 


25.0 

- 2.0 


-10.0 -7.5 -12.0 

3.5 -14.5 3.5 


16.5 

-0.4 


-0.7 
7.0 . 


The lateral accelerations are plotted in Figure 3.4 
for the open loop aircraft, the baseline optimal design, 
and the fixed gain designs. The results of the baseline 
and fixed gain designs are nearly identical. Only in the 
approach flight condition is there any difference in the 
RMS acceleration. For that case, the fixed gain design 
decreases the acceleration slightly. The reason for 
that reduction can be seen in the control surface acti- 
vity show in Figures 3.5 and 3.6. While the different- 
ial flap and separate surface rudder deflections and 
rates of the other four conditions remain nearly the 

same, the deflections and deflection rates increase 

» 

slightly in the approach condition. In the longitudinal 
case an increase is not good because the surfaces are 
working near their RQAS design limits. Here the in- 
creases are not important because the surfaces are below 
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Figure 3.4 Lateral-Directional Control-Rate-Weighting 
Fixed Gain Acceleration Performance 
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Maximum Differential Flap Deflection (deg) Maximum Differential Flap Rate (deg/sec) 



Flight Condition 


Figure 3.5 Lateral-Directional Control-Rate-Weighting 
Fixed Gain Differential Flap Activity 
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their rate and deflection limits. The lateral-direct- 
ional motion does not seem to be affected by flight 
condition as much as the longitudinal motion. 

3.3.3 LONGITUDINAL OUTPUT-WE I GHT I NG FIXED GAINS 


The results of the application of the longitudinal 
output-weighting fixed gain designs were very similar to 
the control-rate-weighting designs. The results can be 
seen in Figures 3.7 through 3.9. The fixed gain that 
was used for this analysis was 


K 


-0.5 

4.5 


0.0 -1.6 -2.5 

0.0 0.5 1.0 


A quick comparison of these gains to the longitudinal 
control-rate-weighting fixed gains of Section 3.3.1 
shows these gains to be quite a bit smaller in magni- 
tude. There are of course only 4 feedback variables 
here as compared to 6 for the control-rate-weighting 
designs. The perturbation velocity feedback terms have 
again been set to zero since the were generally two or 
more orders of magnitude smaller than the other gains. 

The acceleration response of the fixed gain designs 
is nearly the same as the baseline designs for all 
flight conditions. As can be seen in Figure 3.7, only 
the acceleration of the approach flight condition 
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Figure 3.7 Longitudinal Output-Weighting 

Fixed Gain Acceleration Performance 
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Flight Condition 





changes. Here there is a small increase over the base- 
line design. Although the fixed gain design has RMS 
accelerations that are above the RQAS design limit, so 
does the baseline design. The approach flight condition 
was very difficult to design to because the trim flap 
deflection reduced the effectiveness of the flaps as 
direct lift devices. 

The direct lift flap and separate surface elevator 
activity is shown in Figures 3.8 and 3.9. There is a 
wide variation of the surface deflections and rates 
between the baseline designs and the fixed gain designs. 
Both the flap rates and the separate elevator rates were 
within the design limits even with the increases caused 
by the fixed gain designs. The maximum deflections for 
the flaps and separate elevators were not within the 
prescribed limits for several of the flight conditions 
tested. Although this is not good, at the present time 
it is felt that these deflections are tolerable for the 
current aircraft models. If the deflections were li- 
mited, the acceleration would most likely increase. 
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Figure 3.8 Longitudinal Output-Weighting 

Fixed Gain Direct Lift Flap Activity 
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Maximum Elevator Deflection (deg) Maximum Elevator Rate (deg/sec) 



Figure 3.9 Longitudinal Output-Weighting Fixed 

Gain Separate Surface Elevator Activity 


I 
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3.3.4 LATERAL-DIRECTIONAL OUTPUT-WEIGHTING FIXED GAINS 


Once again, the lateral-directional designs were 
much easier to design, and consequently, the fixed gain 
design behave much nicer. The output-weighting lateral- 
directional fixed gain results are given in Figures 3.10 
through 3.12. The fixed gains used are 


K 


2.0 -1.0 -0.7 -1.0 

0.0 1.5 -2.5 1.5 


The zero element in this matrix results from different 
circumstances than the zero elements in the longitudinal 
fixed gain designs. The zero element is the feedback of 
sideslip angle, 6, to the separate surface rudder. This 
would seem to be an important term, and it is. But, in 
this case, the matrix elements from the five baseline 
optimal designs vary around zero giving an average of 
approximately zero. 

The lateral RMS acceleration, given in Figure 3.10, 
shows that the fixed gain designs perform as well or 
better than the baseline optimal designs. The extra 
performance in the approach flight condition once again 
results from increased control surface activity. That 
activity can be seen in Figures 3.11 and 3.12. Although 
there is an increase in differential flap or separate 
surface rudder rates and deflections in some cases, all 


67 


Open Loop 



(2.S/U) SWH Av 


Figure 3.10 Lateral-Directional Output-Weighting 
Fixed Gain Acceleration Performance 
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Flight C 





activity is below the specified RQAS design limits. All 
of the results for the sea-level climb conditions are 
closely matched between the baseline design and the 
fixed gain design. This happens because the fixed gains 
that were chosen happen to be very close to the baseline 
gains for all of the gain matrix positions. 

3.3.5 FIXED GAIN STABILITY 

While the time response of the fixed gain designs 
is very good, there are no assurances that the aircraft 
will be stable under all circumstances. The time simu- 
lations that were performed may not indicate any poten- 
tially dangerous dynamic modes of the aircraft. As a 
further check of the fixed gain response, the eigen- 
values for the closed loop fixed gain aircraft were 
investigated. In addition to indicating the stability 
of the fixed gain controllers for all flight conditions 
and center of gravity locations, the eigenvalues can be 
compared to the eigenvalues of the baseline design to 
give another indication of performance degradation 
caused by the fixed gains. The eigenvalues for the sea- 
level climb flight condition will be discussed below. 
The fixed gain eigenvalues for the remaining flight 
conditions may be found in Appendix B. 
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Maximum Differential Flap Deflection (deg) Maximum Differential Flap Rate (deg/sec) 
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Figure 3.12 Lateral-Directional Output-Weighting Fixed 
Gain Separate Surface Rudder Activity 


71 




The longitudinal control-rate-weighting fixed gain 
eigenvalues are given in Table 3.7 for the sea-level 
climb flight conditon. The eigenvalues for the fixed 
gain system at three center of gravity locations are 
included. While there is some variation of frequencies 
and damping ratios among the center of gravity loca- 
tions, the changes are not large. That is why the time 
simulations for the three center of gravity locations 
all performed nearly identically. The fixed gain phu- 
goid eigenvalues are considerably different than the 
phugoid eigenvalues of the baseline design, but the 
others remain similar. Although the damping on the 
phugoid has been decreased, it is still well damped, 
especially in comparison to the open loop damping. 

Table 3.8 contains the eigenvalue summary for the 
longitudinal output-weighting fixed gain design. Again, 
there is nearly no difference in the eigenvalues when 
the fixed gain designs are implemented at the three 
center of gravity locations. The comparison between the 
baseline eigenvalues and the fixed gain eigenvalues is 
exactly as before. The phugoid increases in frequency 
and decreases in damping when the fixed gains are used. 
The other eigenvalues remain nearly the same. 

The lateral-directional eigenvalues for the sea- 
level climb are given in Table 3.9 for the control-rate- 
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Table 3.7 Longitudinal Control-Rate-Weighting 

Fixed Gain Stability for Sea-Level Climb 


Baseline Eigenvalues: 


Mode 

Frequency 

Damping 

Phugoid 

0.0874 

0.791 

Phugoid 

0.0874 

0.791 

Short Period 

6.02 

0.682 

Short Period 

6.02 

0.682 

Servo 

3.47 

1.0 

Servo 

11.8 

1.0 

Filter 

15.3 

1.0 

Filter 

63.8 

1.0 


i 

I 


Fixed Gain Eigenvalues: 


Foward 

c.g. 

Mid 

Freq. 

Damp. 

Freq. 

0.213 

0.368 

0.214 

0.213 

0.368 

0.214 

6.80 

0.699 

6.44 

6.80 

0.699 

6.44 

2.56 

1.0 

2.88 

11.7 

1.0 

11.7 

15.9 

1.0 

15.7 

63.5 

1.0 

63.7 


g. 

Damp. 

Aft 

Freq. 

c.g. 

Damp. 

0.370 

0.214 

0.372 

0.370 

0.214 

0.372 

0.671 

6.14 

0.643 

0.671 

6.14 

0.643 

1.0 

3.17 

1.0 

1.0 

11.7 

1.0 

1.0 

15.8 

1.0 

1.0 

63.9 

1.0 
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Table 3.8 Longitudinal Output-Weighting Fixed 
Gain Stability for Sea-Level Climb 


Baseline Eigenvalues: 


Mode 

Frequency 

Damping 

Phugoid 

0.0887 

0.886 

Phugoid 

0.0887 

0.886 

Short Period 

11.6 

0.661 

Short Period 

11.6 

0.661 

Servo 

1.94 

1.0 

Servo 

11.3 

1.0 


Fixed Gain Eigenvalues: 


Foward c.g. Mid 


Freq. 

Damp. 

Freq. 

0.216 

0.420 

0.216 

0.216 

0.420 

0.216 

12.1 

0.668 

11.7 

12.1 

0.668 

11.7 

1.77 

1.0 

1.89 

11.3 

1.0 

11.3 


g. Aft c.g. 

Damp. Freq. Damp. 

0.419 0.217 0.418 

0.419 0.217 0.418 

0.658 11.4 0.649 

0.658 11.4 0.649 

1.0 1.99 1.0 

1.0 11.3 1.0 


74 


weighting fixed gain designs, and in Table 3.10 for the 
output-weighting fixed gain designs. The eigenvalues 
for the fixed gain designs and the baseline designs are 
nearly identical for both the control-rate-weighting and 
the output-weighting. Only the servo frequencies for 
the output-weighting show any change. This reinforces 
the findings of the time simulation. There it was noted 
that the fixed gain time response for this flight con- 
ditions was very close to the baseline response. 

3.3.6 GAIN SCHEDULING SUMMARY 

In the preceeding five sections, the time response 
and the stability of fixed gain controllers were inves- 
tigated for five flight conditions and three center of 
gravity locations in both the longitudinal and lateral 
directions. All of the results indicate that the fixed 
gains will perform well enough to prevent the need for 
gain scheduling. There are several things which might 
change these conclusions. If control surface rate and 
deflection limits were strictly enforced, the perform- 
ance of the fixed gains may be degraded beyond accep- 
table levels. There may be an alternative way in which 
to pick gains so that the performance is not degraded by 
the surface limits. Since all of the gains and results 
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Table 3.9 Lateral Directional Control-Rate-Weighting 
Fixed Gain Stability for Sea-Level Climb 



Baseline 

Fixed 

Gain 


Frequency 

Damping 

Frequency 

Damping 

Spiral 

0.594 

1.0 

0.570 

1.0 

Roll 

3.25 

1.0 

3.06 

1.0 

Dutch Roll 

2.75 

0.514 

2.77 

0.507 

Dutch Roll 

2.75 

0.514 

2.77 

0.507 

Servo 

5.47 

0.882 

5.59 

0.871 

Servo 

5.47 

0.882 

5.59 

0.871 

Fi Iter 

12.6 

1.0 

12.4 

1.0 

Filter 

21.3 

1.0 

21.8 

1.0 


Table 3.10 Lateral Directional Output-Weighting 

Fixed Gain Stability for Sea-Level Climb 



Baseline 

Fixed 

Gain 


Frequency 

Damp i ng 

Frequency 

Damping 

Spiral 

0.776 

1.0 

0.683 

1.0 

Roll 

7.39 

1.0 

7.06 

1 . 0 

Dutch Roll 

2.22 

0.578 

2.30 

0.586 

Dutch Roll 

2.22 

0.578 

2.30 

0.586 

Servo 

6.84 

0.874 

7.31 

0.833 

Servo 

6.84 

0.874 

7.31 

0.833 
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depend so directly on the aircraft models, the entire 
topic of gain scheduling and fixed gain designs needs to 
be reinvestigated when the aircraft models are improved. 
If the performance of the fixed gains is close but not 
quite within acceptable levels, it may be possible to 
scheduled a single gain while all of the others remain 
fixed. This is a subject which has not yet been exa- 
mined. 

3 . 4 SUMMARY 

Using the standard optimal regulator solution as a 
basis, the continuous and sampled-data solutions for two 
extended optimal control structures were presented. The 
equations necessary for the solution of these struc- 
tures, control-rate-weighting and output-weighting, were 
presented in a summarized form. Using the design rou- 
tines of ORACLS that have been implemented in the ICAD 
program, the two control structures were applied to a 
Cessna 402B aircraft for five basic flight conditions in 
both the longitudinal and lateral-directional modes. A 
discussion of the complete results for one flight condi- 
tion is included here. Based upon the optimal designs, 
a set of fixed gains were chosen and tested. These 
fixed gains were checked for their feasibilty as a 
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simple solution for the problem of gain scheduling. The 
time response of these fixed gains was compared against 
the response of the baseline designs. With a few ex- 
ceptions, the fixed gains performed well. The excep- 
tions arose from the enforcement of deflection limits. 
The stability of the fixed gains was checked by exam- 
ining the eigenvalues of the closed loop system. The 
fixed gains were all stable and the oscillatory modes 
were all well damped. Based upon these results, the 
fixed gains are a good alternative to gain scheduling. 
This topic will have to reinvestigated as better models 
of the aircraft become available. 


CHAPTER 4. 


SYSTEM MODIFICATIONS 

To continue the investigation of the RQAS , the next 
step is implementation and flight testing of the pro- 
posed control laws. Several hardware related areas need 
to be investigated first. This chapter presents the 
preliminary design of the hardware required to implement 
a RQAS. First, an experimental system, including sen- 
sors and a flight control law computer, is defined. 
Next, an investigation into the possible control surface 
configurations is discussed. The proposed modifications 
to the aircraft control surfaces of the Cessna 402B 
aircraft are then detailed. The control surface actuator 
requirements are also discussed. 

4.1 EXPERIMENTAL SYSTEM 

The basic proposed experimental system, Figure 4.1, 
includes onboard sensors, a digital flight computer, a 
set of analog-to-digital (A/D) and digital-to-analog 
(D/A) converters, an onboard flight engineering station 
and a set of electro-hydraulic actuators, and an onboard 
digital data recording system. No telemetry data record- 
ing equipment is proposed. 
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Figure 4.1 Ride Quality Experimental System 












4.1.1 SENSOR PACKAGE 


A full set of aircraft motion variable sensors is 
required. These include 

1. Accelerometers for vertical and lateral axes 

2. Angle of attack and sideslip angle booms 

3. Vertical gyros for pitch and roll angle 

4. Rate gyros for pitch, roll, and yaw rates 

5. Sensors for temperature and both static and 
dynamic pressure. 

The proposed package includes an angle of attack and 
sideslip angle boom for documentation purposes only. 
This angle data will not be used as control law feedback 
due to the inaccuracy of these types of sensors. In- 
stead, the a and B feedback signals will be based on 
vertical and lateral acceleration data. The temperature 
and pressure sensors are used for, among other things, 
forward velocity determination. Pilot stick positon 
sensors may also be included to aid in the implementa- 
tion of a controller which eliminates problems asso- 
ciated with the RQAS feedback during maneuvers. This is 
an area currently under research. Table 4.1 indicates 
some suggested accuracy and range limits for the sen- 
sors. These values were chosen based upon currently 
available sensors and preliminary estimates of necessary 
data limits. Detailed analysis has yet to be performed. 
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Table 4.1 


Suggested Sensor Requirements 


Variable 

Sensor 

Resolution or 
Accuracy * 

Range 

Ay 

Accelerometer 

0.0020 g's 

+0.5 g's 

A z 

Accelerometer 

0.0024 g's 

+5.0 g's 

e,<t> 

Vertical Gyro 

0 . 5 deg . 

+30 deg. 

P,Q,R 

Rate Gyro 

0.5 deg/sec 

+50 deg/sec 

Temperature Transducer 

2.0 deg. F 

-65 to 120 
deg. F 

Static Pressure Trans. 

0.010 psia 
(25 ft) 

0 to 25,000 
feet 

Dynamic 

Pressure Trans. 

0.005 psia 
(4 knots) 

40 to 150 
knots 


* Whichever value is larger. 


4.1.2 FLIGHT ENGINEERING STATION 

The proposed flight engineering station is composed 
of a mode control panel for controlling the various 
control laws to be implemented and a real time system 
monitoring system. Also, this station will contain the 
digital data recording system. The detailed design and 
selection of equipment for the flight engineering sta- 
tion has yet to be performed. 
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4.1.3 DIGITAL FLIGHT COMPUTER 


The proposed digital flight computer is the ROLM 
1666. This computer is a state-of-the art general pur- 
pose minicomputer designed for airborne applications. It 
is proposed due to its ability to use higher order 
language (FORTRAN) and floating point arithmetic. This 
capability can greatly reduce the time required to de- 
velop and checkout flight coding. The ROLM 1666 has been 
used succesfully in other flight test programs (Refer- 
ence 9). The characteristics of the ROLM 1666 are sum- 
marized in Table 4.2. 
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Table 4.2 ROLM 1666 Characteristics (Reference 9) 


Description: 


Memory : 


General purpose 16-bit minicomputer 
designed to MIL-E-5400 specifications 

32,768 words of 1 psec ferrite core 


Instruction Execution Time: 

(Time in psec, register-to-register operations) 

Fixed Point Floating Point 

(16 bit) (32 bit) 


Add 

Multiply 
Divide 
Load, Store 


1.0 

5.2 to 5.4 

9 . 2 to 9 . 6 

2.0 


1.8 to 4.8 
3.6 to 4.8 
8.0 to 8.8 
4.8 


Support 

Software: Real-Time Disk Operating System (RDOS) 

supports text editors, a FORTRAN 
compiler, and a linking loader, all 
used in constructing executable save 
files. Application programs use the 
Real-Time Operating System (RTOS) for 
flight tests. 


4.2 CONTROL SURFACE MODIFICATIONS 


This section discusses the control surface modifi- 
cations required to implement a ride quality augmenta- 
tion system on a Cessna 402B airplane. The control 
surfaces investigated for the RQAS are a separate sur- 
face elevator and direct lift flaps for longitudinal 
control and a separate surface rudder and differential 
flaps for lateral-directional control. The following 
sections describe the modification options investigated 
and include preliminary structural drawings of the 
recommended options. 


4.2.1 SEPARATE SURFACE ELEVATOR 

To alleviate any feedback of the RQAS control sur- 
face motions to the pilot, an elevator separate from the 
primary elevator is used by the RQAS. The following 
characteristics were assumed as requirements for the 
separate surface elevator (SSE) : 

1. SSE area approximately 20% of the existing 
elevator area. 

2. One SSE surface on each side of the fuselage. 

3. SSE surfaces linked by a single torque tube. 

4. Torque tube driven by a single actuator. 
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Two options which meet these four criteria were 
investigated: 

1. Concentric torque tubes: SSE hinged about the 
primary elevator torque tube (Figure 4.2). 

2. Parallel torque tubes: SSE hinged about an 
axis parallel to and behind the primary 
elevator torque tube (Figure 4.3). 

4. 2. 1.1 CONCENTRIC TORQUE TUBE ARRANGEMENT 

Figure 4.2 illustrates the concentric torque tube 
arrangement. The SSE pivots about the primary elevator 
torque tube, aided by a system of bearings and brackets. 
This arrangement is similar to the separate surface 
elevator system used by Prins et. al. for a separate 
surface stability augmentation system controlled 
Beechcraft Model 99 (Reference 10). 

Although this arrangement offers flexibility in the 
location and size of the separate surface elevator, the 
existing elevator must be modified and the modifications 
are complex. Also, tlje thickness limitations of the 
primary elevator may require a redesign of the primary 
elevator torque tube. 
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Figure 4.2 Separate Surface Elevator: Concentric Torque 
Tube Arrangement 
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Figure 4.3 Separate Surface Elevator: Parallel Torque 
Tube Arranqement 



4 . 2. 1.2 PARALLEL TORQUE TUBE ARRANGEMENT 


Figure 4.3 illustrates the parallel torque tube 
arrangement for the separate surface elevator. The main 
advantage of this arrangement is that the primary eleva- 
tor remains unchanged. The modifications are not complex 
and the hinge moments can be tailored through the posi- 
tioning of the SSE hinge line. However, the size and 
location of the separate surface elevator is restricted 
in this arrangement. 

4. 2. 1.3 SEPARATE SURFACE ELEVATOR PLACEMENT AND SIZING 

Due to the complexity and ultimate cost of the 
concentric torque tube arrangement, the parallel torque 
tube arrangement is considered the most favorable. So as 
not to interfere with the existing elevator, the sepa- 
rate surface elevator is placed in the aft, inboard 
section of the horizontal stabilizer (Figure 4.4). This 
area is unused on the basic 402B. 

Figure 4.5 is an initial structural drawing for 
this arrangement. This drawing is for the SSE section on 
the left side of the fuselage. The total SSE area for 
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Figure 4.4 Chosen Separate Surface Elevator Position 
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Figure 4.5 Preliminary Structural Arrangement of the 
Left Separate Surface Elevator Section 





this configuration is 3.16 square feet or approximately 
21% of the primary elevator surface area. 

The maximum deflection of the separate surface 
elevator is +5 degrees. With the elevator in the maximum 
up position, there will be interference with the rudder 
when it has been deflected approximately 5 degrees in 
either direction. To eliminate this problem, 0.75 inches 
must be trimmed from the bottom of the rudder tab. 
Detailed part drawings for the separate surface elevator 
are included in Appendix C. 

4.2.2 RUDDER MODIFICATIONS 

This section summarizes the design decisions made 
regarding the implementaion of the rudder for use in the 
ride quality augmentation system. The design options 
investigated are: 

1. Separate surface rudder (SSR) . 

2. Ventral fin arrangement. 

3. Twin separate surface rudders located on the 
horizontal stabilizer. 

4. Use of the entire rudder 

(no structural modification). 

The following sections point out the characteris- 
tics of each option and discuss the final decision. 
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4. 2. 2.1 SEPARATE SURFACE RUDDER 


The separate surface rudder was investigated as an 
option which, like the separate surface elevator option, 
would not feedback the ride quality augmentaion system 
motions to the pilot. It was found in the RQAS designs 
that a separate surface rudder with an area equal to 
approximately 33% of the area of the existing rudder is 
needed for adequate control power for the RQAS. 

Two locations of the SSR were investigated. The 
first placed the SSR at the top of the primary elevator; 
the second placed the SSR at the bottom of the primary 
rudder . 

Figure 4.6 depicts the placement of the SSR at the 
top of the primary rudder. This arrangement is compar- 
able to the arrangement used in the SSSA airplane of 
Reference 10. An actuator in the vertical tail drives 
the SSR through an arm fixed to the SSR. This arrange- 
ment would be relatively easy to actuate and implement 
but would require a redesign of the mass balance of the 
primary rudder. 

Figure 4.7 depicts the placement of the SSR at the 
bottom of the primary rudder. This arrangement requires 
no modification to the primary rudder mass balance but 
would require modifications to the rudder trim system. 
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Also, mounting would be difficult since the primary 
rudder torque tube translates as it rotates. 

In both location options, a large portion of the 
primary rudder would be lost to the separate surface 
rudder. This could cause problems with directional con- 
trol and may require the addition of some mechanism to 
allow the use of the SSR in conjunction with the primary 
rudder during takeoff and landing. 

4. 2. 2. 2 VENTRAL FIN ARRANGEMENT 


Another configuration which would provide a surface 
for directional control which is separate from the pri- 
mary rudder is a ventral fin arrangement. Figure 4.8 
illustrates this type of configuration. The low aspect 
ratio of the fin results from a compromise in required 
fin area (approximately 33% of the rudder area) and 
take-off rotation angle. The aft fuselage angle is re- 
duced by 2 degrees. 

The main advantage of this arrangement is the fact 
that it requires no modification to the existing rudder. 
The disadvantages are that the fin would be cumbersome 
to deflect in flight and the low aspect ratio would 
provide uncertain control power. Also, as noted above, 
the ventral fin reduces the take-off rotation angle. 
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4. 2. 2. 3 TWIN SEPARATE SURFACE RUDDERS 


The twin separate surface rudder arrangement is 
illustrated in Figure 4.9. This arrangement is again 
designed to provide directional control separate from 
the primary rudder. The twin "rudders" act like variable 
incidence vertical stabilizers. 

This arrangement again requires no modification to 
the existing rudder. However, it does require the stiff- 
ening of the horizontal stabilizer structure and may 
cause interference with the airflow over the existing 
elevator making the elevator control power unpredict- 
able. 

4. 2. 2. 4 USE THE ENTIRE EXISTING RUDDER 

After examining the options which provide direct- 
ional control for the RQAS separate from the primary 
rudder system, it is seen that each has major drawbacks. 
For this reason, it was decided to devote the entire 
existing rudder to the RQAS. The fact that this option 
requires no modifications to the rudder is the prime 
factor in this decision. The negative aspect is that the 
pilot will now have to deal with the feedback from the 
RQAS. 
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Since this is a research project, it is felt that 
the feedback will not be a major problem since the pilot 
will be warned of the problem and much of the flight 
testing will be performed with the pilots feet on the 
floor. Also, the RQAS rudder deflections will be limited 
to +5 degrees. If the RQAS were to be incorporated into 
the design of a new aircraft, the separate surface 
rudder option could be easily implemented and would be 
recommended. 

4.2.3 FLAP MODIFICATIONS 

As stated earlier, the flap must perform two oper- 
ations. The entire flap is driven in the longitudinal 
mode. In addition, the outboard halves of the flap are 
driven in the lateral-directional mode. Also, the flap 
must be modified to deflect up 15 degrees and down 45 
degrees. The present flap sytem on the Cessna 402B is a 
split flap system and only deflects down. Two basic 
arrangements were investigated: 

1. Modify the wing to include plain flaps. 

2. Use a trailing flap system (external airfoil). 
The following sections describe each of these options in 
more detail. 
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4. 2. 3.1 PLAIN FLAP MODIFICATION 


The plain flap modification option requires rework- 
ing the trailing edge of the wing, aft of the rear spar. 
The piano hinge on the lower spar cap is used to attach 
the flap to the spar (Figure 4.10). No major structural 
modification is required and the inboard and outboard 
sections of the flap can easily be actuated separately. 
Actuating the outboard section separately is a require- 
ment due to the dual role of this section. The disad- 
vantage of this modification is the fact that the wing 
locker must be modified to allow the upward deflection 
of the flap. This leaves a section of the engine nacelle 
extending beyond the trailing edge of the flap (Figure 
4.11) and makes it hard to predict the control power of 
the flap. Also, the actuation mechanism must be rede- 
signed. 
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4. 2. 3. 2 TRAILING FLAP SYSTEM 


As a simple alternative to the plain flap modifica- 
tion, a trailing flap (auxiliary airfoil) could be 
mounted below and aft of the wing trailing edge (Figures 
4.11 and 4.12). However, several problems arise from 
this arrangement. First, additional structural members 
are required for supporting the flap. Also, this ar- 
rangement may cause severe changes in the center of 
gravity and stability characteristics of the airplane. 

4.2.3. 3 FLAP DESIGN CHOICE 

For this project, the wing mounted plain flap con- 
figuration was chosen. The basic configuration is il- 
lustrated in Figures 4.13 and 4.14. These figures illus- 
trate the inboard and outboard sections of the proposed 
flap, respectively. This arrangement allows the outboard 
flap sections to be deflected dif ferentialy and the 
inboard flap sections to be deflected symetr ically . 

As noted earlier, this configuration requires some 
modification to the wing locker portion of the engine 
nacelle. Figure 4.15 illustrates the required modifica- 
tions. Detailed part drawings for the flap modifications 
are included in Appendix C. 
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Figure 4.12 Trailing Flap System (External Airfoil) 
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Figure 4.13 Inboard Flap Preliminary Structural Drawing 
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4.2.4 SUMMARY OF PROPOSED RQAS CONTROL SURFACE 
MODIFICATIONS 


The proposed control surface modifications to the 
Cessna 402B for the purposes of the ride quality aug- 
mentation system are summarized as follows: 

Longitudinal Mode 


1. Separate surface elevator located in the aft, 
inboard section of the horizontal stabilizer 
with a deflection range of +5 degrees. 

2. Plain flap replacing existing split flap with 

a dif ferent ialy deflecting outboard section and 
a symetrically deflecting inboard section. Each 
section can deflect +15 to -45 degrees. 


Lateral-Directional Mode 


1. Outboard section of the plain flap described 
above . 

2. Use the entire existing rudder (limiting the 
RQAS range of deflections to +5 degrees). 

It should be remembered that the modifications 
proposed were motivated by providing the capability to 
evaluate the RQAS concept with minimum modification to 
the vehicle. The configuration could be significantly 
different in new aircraft designs. 


If additional modes are required, e.g. engine out 
control, these limits may be changed. 
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4.3 CONTROL SURFACE ACTUATION 


This section summarizes the preliminary actuator 
sizing for the Cessna 402B ride quality augmentation 
system. The sizing involves calculation of the maximum 
hinge moments for each control surface and the corres- 
ponding maximum actuator loads, stroke, and speed. A 
hydraulic actuation system is proposed with actuators 
controlled by electric servo valves. The components 
needed for this system are a hydraulic pump to be mount- 
ed on one of the engines, an accumulator and a set of 
electro-hydraulic actuators. The proposed system will 
operate at hydraulic pressures from 1,000 to 3,000 
pounds per square inch. 

4.3.1 RQAS CONTROL SURFACE GEOMETRY 

The geometry of the control surfaces proposed for 
the ride quality augmentation system in Section 4.2 is 
summarized in Table 4.3. This information is needed in 
the hinge moment calculations. The separate surface 
elevator information is given for the sum of the two 
surfaces since they will be driven by a single actuator. 
The flap information is for the Siam of the inboard 
sections which will also be driven by a single actuator. 
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The differential flap information is for a single sur- 
face since the left and right sections will be driven 
separately. The rudder requires a single actuator. 


Table 4.3 

RQAS Control 

Surface Geometry 


Surface 

Span 

S(ft) 

Mean Chord 
c ( f t ) 

Area 

S(ft 2 

SS Elevator 

2.66 

1.19 

3.16 

Flap( Inboard) 

7.43 

1.53 

11.36 

Differential Flap 

4.95 

1.50 

7.41 

Rudder 

6.67 

2.66 

17.77 


4.3.2 HINGE MOMENT CALCULATIONS 

This section summarizes the hinge moment calcula- 
tions used in the actuator sizing. The hinge moment 
coefficients were obtained from Reference 3 and were 
verified by DATCOM methods (Reference 11). The hinge 
moments were calculated using the equation: 

HM = Ch q b c 2 (4.1) 

where HM is the hinge moment, Ch is the hinge moment 
coefficient for the given control surface, q is the 
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design dynamic pressure (for V max at sea level), and b 
and c are the control surface span and mean chord re- 
spectively (Table 4.3). The hinge moment calculations 
are summarized in Table 4.4. 

Table 4.4 Hinge Moment Summary 


Surface 

Hinge Moment 
Coefficient 

Hinge Moment 
Equation 

SS Elevator 

0.0219 

178 

Flap( inboard) 

0.0450 

1,691 

Differential Flap 

0.0450 

1,083 

Rudder 

0.0145 

1,478 


4.3.3 ACTUATOR REQUIREMENTS 

The characteristics of the required actuators are 
summarized in Table 4.5. The calculations were made 
assuming a 4 inch moment arm and a linear actuator for 
each control surface. The specifications given in Table 
4.5 can be modified easily by varying this moment arm. 
The maximum load calculations were made using the fol- 
lowing equation: 

Max Load = 1.4 x HM( in-lb) /moment arm(in) » (4.2) 
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where the 1.4 is a factor of safety. The speed and 
stroke are based on the maximum control surface rate and 
deflection respectively. The speed is calculated with 
the equation: 

• 

Speed = xh, (4.3) 
• 

where 6 max the maximum control surface rate and h is 
the moment arm. The stroke is based on the geometry of 
the mechanism and the deflection range. 
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Table 4.5 Actuator Requirements 


Control Surface Deflection and Rate Limits 


Surface 

Deflection 
Range (deg) 

Maximum 

Rate(deq/sec) 

SS Elevator 

+5 

50 

Flap( inboard) 

+15 to -45 

120 

Differential Flap 

+15 to -45 

120 

Rudder 1 

+ 32 

50 


Actuator Requirements 


Surface 

Max Load 
(lbs) 

Speed 
( in/sec) 

Stroke 
( in) 

SS Elevator 

65 

3.50 

0.75 

Flap( inboard) 

750 

8.50 

4.25 

Differential Flap 

380 

8.50 

4.25 

Rudder 

520 

3.50 

4.50 

1 The deflection 
uses a deflection 

is for the 
range of +5 

standard rudder, 
degrees. 

The RQAS 
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CHAPTER 5. 


SUMMARY, CONCLUSIONS, AND RECOMMENDED RESEARCH 
5.1 SUMMARY 

The primary goals of this phase of work have been 
the generation of detailed optimal designs for the long- 
itudinal and lateral-directional RQAS. Another goal was 
the preliminary investigation of the necessary hardware 
modifications to the Cessna 402B to implement separate 
surface controllers and direct lift flaps. These goals 
have been achieved along with the groundwork for the 
implementation of the RQAS controller on an experimental 
aircraft. This groundwork included the preliminary 
definition of the experimental system and the deter- 
mination of the necessary requirements for surface 
actuation. 

The results of this (and previous) work show that 
substantial reductions can be made in the RMS values of 
the vertical and lateral accelerations caused by atmos- 

t 

pheric turbulence by implementing an active control 
system. It was shown in Section 3.3 that the implemen- 
tation of fixed gain designs, based on either the con- 
trol-rate-weighting or output-weighting optimal design 
techniques, give good results across the basic flight 
envelope. 


113 


5 . 2 RECOMMENDED RESEARCH 


The next steps in the RQAS project should involve 
all of the factors necessary to implement and experi- 
mental flight test the system. The following is a list 
of research topics which follow from the present pro- 
ject. Some are directly related to the implementation 
of the experimental system and others relate to fur- 
thering the basic design concept. 

1. Continue with the hardware modification plans 
including the detailed structural design of the flaps 
and separate elevator. This includes a detailed 
structural analysis along with flutter analysis. Of 
special concern is the bending moments caused by the 
new, more active, control surfaces. 

2. Proceed with detailed study of the experimental 
system and hardware^ requirements for implementation 
of the controller on the Cessna 402B. This includes 
investigation of sensor accuracy and placement, digi- 
tal computer requirements, and data recording equip- 
ment . 

3. Continue with the analytic designs using the 
RQAS concept. This might include investigations of 
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other control structures such as Limited State Feed- 
back. 

4. Investigate the response and interaction of the 
individually designed longitudinal and lateral-direc- 
tional RQAS controllers on fully coupled linear mo- 
dels. 

5. Investigate methods to prevent interference with 
the pilot induced accelerations. This might be a- 
chieved by including washout filters in the control 
circuits or by directly measuring the pilot commands 
and implementing a trim map to account for the ac- 
celerations requested. 

5. Study the effects of unsteady aerodynamics on 
the RQAS performance. The effects include the lag of 
control response when the control surfaces are de- 
flected at high rates, the effects of the unsteady 
flap wake on the tail surfaces and controls, and the 
structural interactions caused by the control acti- 
vity. 

7. Develop a piloted simulation using the C-402 B 
model operational at NASA Langley. This simulation 
would be used to evaluate RQAS designs, conduct fail- 
ure mode analysis, and evaluate candidate techniques 
for RQAS operation during intentional maneuvering. 
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8. Perform parameter identification flight test of 
the modified vehicle. If major model variations 
occur, refine the RQAS algorithms. 

9. Perform a flight test program of a digital RQAS. 
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APPENDIX A 


EQUATIONS OF MOTION WRITTEN IN STATE-SPACE FORM 

This appendix gives the conversion of the standard 
aircraft small perturbation equations of motion found in 
Reference 4 to state-space form. The equations are writ- 
ten in the stability axis system and are modified 
slightly to account for the change in the acceleration 
vector due to its rotation through the angle of attack. 
The resulting equations are consistent with the models 
supplied by NASA LaRC for the Cessna 402B aircraft 
(Reference 3 ) . 

The equations are separated into the longitudinal 
and lateral-directional modes and the modified dimen- 
sional stability derivatives are defined. 


A. 1 


Assuming steady-state straight and level trimmed 


flight: 

Vi = «j>i = w x = 0 ; 

Pi = Ql = Rl = 0 ; 

• • • 

♦l ■ 01 * $1 * 0 ; 

the equations of motion as derived in Reference 4 can be 
written as: 

mu = -mgycosYi + f A + f T 

A A 

m(v + U^r) = mg<J>cos6i + f A ^ + f^ 

m(w - Ujq) = -mgysinyi + f A + fT 

z z 

• • 

*xxP " I xz r = 1-A + 

I yyCJ = 

• • 

Izz r ~ *xzP = n A + n T 


p = 4> — ipsin9i 
q = 0 


r = ipcosQi 


The small perturbation equations of motion can then 
be written in the dimensional stability derivatives as 
follows. 


Longitudinal : 

• • 

w - Uiq * -gysinyi + Z u u + Z^a + Z a a + Z q q + 

+ + Z§^6f 

u = -gycosYi + X u u + X a a + X$ 6 se + X$ 6f 

S 6 I 

( A. 1 ) 

* . • 

q = M u u + M a a + M a a + Mqq + Mg 6 se + Mg Sf 

S 6 t 

e = q 


Lateral-Directional : 

v + Uir = g<J>cos0i + Y p p + Yg0 + Y r r + Yg df 6 df + 

• • 

p - Air = LgB + L p p + L r r + Lg df 6 d f + Lg sr> 6 sr 


(A. 2) 


r - Bip = N p p + NgS + N r r + Ng df 6 df + 


N« sr «sr 


A. 3 


♦ = p + tanQir 



where the 4> equation has been neglected and 

Ai = Ixz/^xx and = Ixz/^zz 

The dimensional stability derivatives as defined in 

Reference 4 are presented in Tables A.l and A. 2. 

Introducing the equations 

* • 
w = U]_a 

v = UiB 


Y = 0 - a 


into Equations (A.l) and (A. 2) and rearranging results 
in the following state-space representation of the small 
perturbation equations of motion where the modified 
dimensional stability derivatives are defined in Tables 
A. 3 and A. 4. 


Longitudinal : 


a 

u 

q 

6 


Za 

*a 

Mi 


0 


Zi Z4 z& 


a 


Z *se z *f 

Xi 0 x$ 


u 

+ 

XJ se x if 

Mi M$ 


q 


“ise MJ f 

0 0 0 


0 


1 

o 

o 

1 


$se 

«f 


(A. 3) 


A. 4 


Lateral-Directional: 




A. 5 


Table A. 1 Longitudinal Dimensional Stability 
Derivatives (from Reference 4) 


X = 
u 


*T 


-’l S < C D + V 

u (sec ) 


mU, 


q 1 s(c T + 2c x ) 


mU 


x x i _i 

u — (sec ) 


q . ScC 
M 1 m 


M = — r- 
a I 


yy 


-2 

(sec ) 


-’l S(C D - C L> 


a m 

-il SC D 


^ (ft sec 


6 _ /r -2 
E (ft sec or 


6 e ra 


r -2. -lv 

ft sec deg ) 


V‘ C L + 2 C L> 

z - rtj (seC * 

u mU- 


Hj, = i 

a yy 


q.Sc 2 C 
M 1 m* 


M- = 


a 21 U. 

yy i 


(sec Z ) 


(sec S 


il S ( C L + C D ) 

Z - (ft sec '■) 


- 2 . 


m 


q Sc z C 

m 


M = 


q 21 U. 

4 yy i 


(sec 


4 i sc l- c 


Z 'a " 2mU. 


V C L C 


Z q 2mU 


(ft sec 1 ) 


(ft sec ■*■) 


q. ScC 
H 1 m. 


M r 


/ “2 
E (sec or 

1 -2 , -I, 

yy sec deg ) 


4 i sc l, 


_2 

J E (ft sec or 


m 


_ “2, ~U 

ft sec deg ) 


M = 
u 


q.Sc(C + 2C ) 
n l m m 1 , 

(ft sec ) 

I U. 

yy 1 


"is 


«i Sc(c «. T + 2 \ ) 

„ Li i- (ft* 1 sec -2 ) 

\ 
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I cr | l cr | I cr 


Table A. 2 Lateral-Directional Dimensional Stability 
Derivatives (from Reference 4) 


<i sbc t, 


1 (ft sec' 2 ) 


L SbC 

r— 2. (ft- sec ■*•) 

2mU^ 


. SbC 

1 y r -1 

Y r = “Suf (ft S6C } 


q, SC _ 2 

y S. (ft sec or 

Y = - -2 -1 

6, m ft sec deg ) 

A 


q SC _ 2 

y 6 (ft sec or 

Y = 5 — -2 -1 

6 m ft sec deg ) 

R 


q i SbC t 

x -2 
Lg l (sec Z ) 

XX 


q l sb2c t 

L p = — uf (sec > 

r XX 1 


q l sb2c * 


l r ■ 21 u/ (serl) 

XX 1 


A , -2 

(sec or 

xx -2 , -lv 

sec deg ) 


q i Sbc n, 


R , -2 

— (sec or 


xx -2 , -1,. 

sec deg ) 


q, SbC 

1 n R -2 

Ng = — j (sec ) 

zz 


q, SbC 
1 n T 

N = — (sec 2 ) 

T S 'zz 


q, Sb 2 C 
1 n . 

N P = ~n uT 

r ZZ 1 


2- (sec ■*■) 


q l Sb2C n . 

»,-rr (sec } 

zz 1 


q. SbC 

1 n A 

°A -2 

n = — (sec or 

o X — 2* <— 1 

" 22 sec deg ) 


q, SbC 

1 n* 

°R -2 

N a — (sec or 

R ZZ cpp Z dp 


-2 . - 1 , 
sec deg ) 
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Table A. 3 Modified Longitudinal Dimensional Stability 
Derivatives (from Reference 7) 


Zi = (gsinYi+Z a )/(Ui-Za) 
Zu “ Z u /(U].-Z a ) 

Z^ = (Ui+Z q )/(Ui-Za) 

Z§ = -gsinYi/(Ui-Za) 

Z k = z 5 /(Ui-zi) 

se se 

z$ f = z 6f /(ux-z;) 

Xi = gcosYi + x a 
xi = x u 


x$ = x$ 

se u se 
x i f = X 6f 

Mi = M^i + M a 
Mi = M^Zi + M u 
Mq = M^i + M q 

ms = m;z§ 

M4 = M^Z* + m 6 
se se u se 

Mi f = MaZ4 f + M 5f 


Xi = -gcosYl 


Table A. 4 Modified Lateral-Directional Dimensional 

Stability Derivatives (from Reference 7) 


Y& = Y B /Ui 
Y p = Yp/Ui 
Yf = (Y r /Ui)-1 
Y| = gcos8i/Ui 

Y^df = Y ^df /Ul 

Y *sr = Y ^sr /U l 

L£ = (A^Ng+Lg) / ( l-A^B^) 
Lp = (AxNp+LpJ/d-AiB!) 
L r = (AiNr+LrJ/d-AiBi) 


L^ df = (A 1 N 5df+ L 6df )/(l-A 1 B 1 ) 

L^ sr - (AiN 6sr+ L 6sr )/(l-A 1 3 1 ) 

Nfc = (B 1 Lg+Ng)/(l-AiB 1 ) 

Np = (B 1 L p +Np)/(l-A 1 B 1 ) 

Nf = (BiL r +N r )/(l-AiBi) 

= (BiL 6df +N 6df )/(l-Ai3 1 ) 

Ni sr = (Bl L 5 sr +N 6 sr ^( 1 - A l 3 l) 



APPENDIX B. 

MATHEMATICAL MODELS AND OPTIMAL DESIGN RESULTS 

This appendix presents the basic math models of the 
Cessna 402B. These are perturbation models that have 
been linearized about five basic flight conditions as 
described previously. The models are presented here in 
terms of four basic matrices which satisfy the fol- 
lowing linearized equations 
x = Ax + Bu 
y = Cx + Du 

The models for 3 center of gravity locations in the 
longitudinal direction are presented as are the basic 
lateral-directional models. 

Following the tables of the aircraft design models 
are the tables of design results. These tables give the 
complete results of the baseline designs for all flight 
conditions. The results for both control-rate-weighting 
and output-weighting are given. The tables, beginning 
at Table B.21, give the basic time response of each 
design along with the closed loop eigenvalues, the op- 
timal gain matrix, and the diagonal elements of the 
optimal weighting matrices. 

The stability of the fixed gain models is shown in 
Tables B. 41 through B.60. The eigenvalues of the base- 
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line designs are compared to the eigenvalues for the 
fixed gain designs in all of the flight conditions and 
at all of the center of gravity locations previously 
investigated. Both the control-rate-weighting and out- 
put-weighting designs are checked for stability. 
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Table B.l Cessna 402B Longitudinal Mathematical Model 
For Sea Level Take-off (CG = 0.25) 


X 

= Ax 

+ Bu 





• 

a 


i -1.1730 

| 

-0.0017 

0.9133 

-0.0249 

r 1 

: a ! 

i 

• 

u 


j 9.6588 

-0.0278 

0.0000 

-31.7843 

j 

u ; 

i : i 

• 

q 


i -5.4978 

i 

0.0007 

-7.5327 

0.0784 1 | 

j ; 

q 

• 

0 


0.0000 

0.0000 

1.0000 

0.0000 


e i 





— 







-0.0361 

-0.2262 






0.0000 

-4.5082 


^se 




+ 

-3.7650 

1.1163 


_ 6f _ 





0.0000 

0.0000 



y = Cx + Du 


Hi 

1 


-220.3845 

-0.3190 

-15.9281 

0 . 1 3 7 7~ 


— — 









a 


a 


1.0000 

0.0000 

0.0000 

0.0000 


u 


u 

— 

0.0000 

1.0000 

0.0000 

0.0000 


q 


q 


0.0000 

0.0000 

1.0000 

0.0000 


e 


e 


0.0000 

0.0000 

0.0000 

1.0000 







— 

* 










”-6.6301 

-41.5912 








0.0000 

0.0000 


i — 

! 

-i 








$se 




+ 

; o.oooo 

0.0000 


s f 






; 0.0000 

0.0000 



J 


! 0.0000 0.0000 
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Table B.2 Cessna 402B Longitudinal Mathematical Model 
For Sea Level Climb (CG = 0.25) 


x = Ax + Bu 


• 





- 

- 

a 


-1.3325 

-0.0014 

0.9189 

-0.0120 

i 

• 

u 


12.7885 

-0.0228 

0.0000 

-32.0688 

i 

j 

• 

q 


-6.4781 

0.0023 

-8.1525 

0.0406 


• 

6 

1 

0.0000 

0.0000 

1.0000 

0.0000 

- 


-0.0389 

0.0000 

-4.6678 

0.0000 


-0.2595 

-5.9354 

1.4135 

0.0000 



y = Cx + Du 


— 1 





— 


— - 


a z 


-283.5580 

-0.2857 - 

-17.1010 

0.0780 

i 

I 

a 


a 


1.0000 

0.0000 

0.0000 

0.0000 


u 


u 


0.0000 

1.0000 

0.0000 

0.0000 


q 


q 


0.0000 

0.0000 

1.0000 

0.0000 


0 


e 


0.0000 

0.0000 

0.0000 

1.0000_ 




— 


“ 

* 










-8.2010 

-54 . 7130*1 








0.0000 

0.0000 


r 

1 








$se 




+ 

0.0000 

0.0000 


6 f 






0.0000 

0.0000 








0.0000 

0.0000 
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Table B . 3 Cessna 402B Longitudinal Mathematical Model 
For Climb at 5000 Feet (CG = 0.25) 


x = Ax + Bu 


-1.2413 

12.5650 

-7.1464 

0.0000 


■ 0.0012 

■ 0.0212 

0.0016 

0.0000 


y = Cx + Du 


1 — — 

a 2 

r 

a 


u 

i 


! 

q 1 

: 

f 

CD j 



-284.6474 
1.0000 
0.0000 
0.0000 
_ 0.0000 


0.9304 

0.0000 

■7.5786 

1.0000 

-0.0361 

0.0000 

•4.6749 

0.0000 


- 0.01051 


-32.0750 u 
0.0330 i q 


o.ooooj l^ej 


-0.2417 I 

-5.9354 | 5 se 

1.2914 6f 
L 1 J 

0.0000 


-0.2662 

0.0000 

1.0000 

0.0000 

0.0000 


-15.8299 

0.0000 

0.0000 

1.0000 

0.0000 

-8.2155 
0.0000 
0.0000 
0.0000 
_ 0.0000 


0.0636 

0.0000 

0.0000 

0.0000 

1.0000 

-54.9407 

0.0000 

0.0000 

0.0000 

0.0000 
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Table B. 4 Cessna 402B Longitudinal Mathematical Model 
For Cruise at 20,000 Feet (CG = 0.25) 


x = Ax + Bu 


• 

: 

— 



“t 1 

a 


-1.2343 

-0.0005 

0.9583 

o . oooo : 

i 

• 

u 


18.3366 

-0.0178 

0.0000 

-32.1733 


• 

q 


-11.7142 

0.0008 

-7.1964 

0.0000 


• 

e 


0.0000 

0.0000 

1.0000 

0.0000 



■0.0343 

0.0000 

•6.9711 

0.0000 


-0.2423 
-9.2750! [fi S( 

1.5455! 6f 

! L 

0.00001 


-441.7867 

1.0000 

0.0000 

0.0000 

0.0000 


■0.1736 

0.0000 

1.0000 

0.0000 

0.0000 


-14.9157 

0.0000 

0.0000 

1.0000 

0.0000 


0.0000 

0.0000 

0.0000 

0.0000 

1 . 0000 _ 


-12.2596 -86.7313 

0.0000 0.0000 

0.0000 0.0000 

0.0000 0.0000 

0.0000 0.0000 
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Table B.5 Cessna 402B Longitudinal Mathematical Model 
For Sea Level Approach (CG = 0.25) 


x = Ax + Bu 


• 


— 



- 

r 

- n 


a 


-1.0456 

-0.0024 

0.9192 

0.0102 


a 


• 

u 


21.0420 

-0.0525 

0.0000 

-32.1253 

t 

u 


• 

q 


-3.7424 

0.0049 

-6.2240 

-0.0264 1 

q 


• 

e 

„ __ i 


0.0000 

0.0000 

1.0000 

0.0000 

-* L 

e 




t 

1 

-0.0297 

"1 

-0.0570 






0.0000 

-3.1464 


6 S e 



+ 




i 





-2.7075 

1.0004 






0.0000 

0.0000 




J 


y = Cx + Du 


a z 


-165.8660 

-0.3822 

-12.9420 

-o.osir 


— 

a 


1.0000 

0.0000 

0.0000 

0.0000 


cx 








u 

u 

= 

0.0000 

1.0000 

0.0000 

0.0000 










q 

q 

j 


0.0000 

0.0000 

1.0000 

0.0000 


i a 

H 


0.0000 

0.0000 

0.0000 

1.0000 

t J 


-4.7561 

-9.1360 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000_ 


$se 

5 f 


+ 


Table B.6 Cessna 402B Longitudinal Mathematical Model 
For Sea Level Take-off (CG = 0.34) 


x « Ax + Bu 


■1.1739 

9.6588 

-1.8837 

0.0000 


•0.0017 

-0.0278 

0.0004 

0.0000 


0.9157 

0.0000 

-7.1388 

1.0000 

-0.0360 

0.0000 

-3.6665 

0.0000 


-0.0249 a 
-31.8268 u 
0.0743 q 

i 

I 

o.oooo I e 


-0.2264 

-4.5082 

1.8218 

0.0000 


-215.8356 

-0.3126 

-15.4996 

3.8543 : 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000_ 


a 

' 

i 

! u 


-6.6264 

0.0000 

0.0000 

0.0000 


!_ 0.0000 


-41.6264 

0.0000 

0.0000 

0.0000 

0.0000 
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Table B. 7 Cessna 402B Longitudinal Mathematical Model 
For Sea Level Climb (CG = 0.34) 


X 

= Ax ■» 

Bu 






r- 

• 


— 



- 


r— 


a 


- 1.3335 

- 0.0014 

0.9210 

- 0.0120 




• 

u 


12.7885 

- 0.0228 

0.0000 

- 32.0688 


u 


• 

q 


- 1.7499 

0.0021 

- 7.7263 

0.0385 


q 


• 

0 


0.0000 

0.0000 

1.0000 

0.0000 


e 


“ — J 

i 

- 



-i 

— 

i 

k 




i 

i — 







1 

|- 0 . 0389 

- 0.2597 







| 0.0000 

- 5.9354 


o* 

w 

!? 1 




+ 

1 

1 

- 4.5455 

2.3424 








0.0000 

0.0000 




y = Cx + Du 


a z 


-281.1698 

-0.2952 - 

-16.6572 

2.8614 - 

' — 









a 


a 


1.0000 

0.0000 

0.0000 

0.0000 


u 


u 

= 

0.0000 

1.0000 

0.0000 

0.0000 




q 


0.0000 

0.0000 

1.0000 

0.0000 

i 

i 

q 








i 

e 


L e j l 

0.0000 

0.0000 

* 

0.0000 

1.0000 

1 

— 





-8.2063 

-54.7580 







0.0000 

0.0000 

! 

r 




+ 

0.0000 

0.0000 

I 

5 se 




0.0000 

0.0000 

1 

L«f 

j 




0.0000 

0 . 0000 _ 
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Table B.8 Cessna 402B Longitudinal Mathematical Model 
For Climb at 5000 Feet (CG = 0.34) 


x = Ax + Bu 

-1.2421 -0.0012 0.9322 -0.0105 

12.5650 -0.0212 0.0000 -32.0808 

-2.3755 0.0014 -7.1815 0.0313 

0.0000 0.0000 1.0000 0.0000 

t — “i 

-0.0361 -0.2419 

0.0000 -5.9355 6 se 

+ 

-4.5520 2.2285 6f 

0.0000 0.0000 

— i 

y = Cx + Du 





-282.3740 

-0.2728 

1.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 


-15.4134 

2.8603 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

1.0000_ 



f-8 .2205 

I 0.0000 

i 

+ j 0.0000 
0.0000 
0.0000 


-54.9926 

0.0000 

0.0000 

0.0000 

0.0000 



B.10 


Table B.9 Cessna 402B Longitudinal Mathematical Model 
For Cruise at 20,000 Feet (CG = 0.34) 


x = Ax + Bu 


• 


— 



— 


— . — 

a 


-1.2348 

-0.0005 

0.9594 

0.0000 


a 

• 

u 


18.3366 

-0.0178 

0.0000 

-32.1740 


u 

• 

q 


-4.1742 

0.0007 

-6.8172 

0.0000 


q 

• 

e 


0.0000 

0.0000 

1.0000 

0.0000 i 


0 

- — 


- 



— i I 

- -i 


i 

| y = Cx + Du 


j-0.0343 

0.0000 

-6.7860 

0.0000 


-0.2424 

-9.2750 

3.0365 

0.0000 


5 se] 
5 f ! 


a z 

a 

u 

q 

e 


-441.9510 

-0.1790 

-14.5313 

1.5492 

c - - 






j 

a 


1.0000 

0.0000 

0.0000 

0.0000 

i 

! 

u 


0.0000 

1.0000 

0.0000 

0.0000 

; 







i q 


0.0000 

0.0000 

1.0000 

0.0000 

! 

e 


0.0000 

0.0000 

0.0000 

1.0000 





L 


— 





1 

! 

j-12 .2693 

-86.7581 





! 

! 0.0000 

0.0000 


r 




! 



!$se 


+ ’ 

0.0000 

0.0000 


5 f 




0.0000 

0.0000 






! 

1 0.0000 

0.0000 
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Table B.10 Cessna 402B Longitudinal Mathematical Model 
For Sea Level Approach (CG = 0.34) 


X 

= Ax 4 

Bu 





— - 







— — 1 

a 


- 1.0464 

- 0.0024 

0.9213 

0.0102 ; 

i 

a | 

i 

• 

u 


21.0420 

- 0.0525 

0.0000 

i 

- 32.1299 

u 

• 

q 


- 0.9632 

0.0046 

- 5.8986 

- 0.0250 

q 

• 

e 


0.0000 

0.0000 

1.0000 

0.0000 

0 : 






— * 

— i 





— 







- 0.0297 

- 0.0570 ; 

! 






0.0000 

- 3.1464 j 

$se 




+ 

- 2.6365 

1.1559 j 

J 

«f _ 





0.0000 

0.0000 ! 
J 



y = Cx + Du 


HaJ ) 1 

-167.6825 

-0.3846 - 

■12.6114 

-0.5161 

1 — 1 








! a 


a 


1.0000 

0.0000 

0.0000 

0.0000 

1 

i 

! u 


u 

= 

0.0000 

1.0000 

0.0000 

0.0000 


q 


q 


0.0000 

0.0000 

1.0000 

0.0000 


e 


6 


0.0000 

0.0000 

0.0000 

1 . 0000 J 






L 








■ 

-4.7593 

-9.1341 






0.0000 

0.0000 

r~ 

~1 




l 



: *se 



+ 

0.0000 

0.0000 


1 

1 








«f 




0.0000 

0.0000 


l- 

. — 




0.0000 

o . oooo _ 




B. 12 


Table B.ll Cessna 402B Longitudinal Mathematical Model 
For Sea Level Take-off (CG = 0.14) 


X 

= Ax - 

► Bu 






• 

a 

I 

- 1.1719 

- 0.0017 

0.9106 

- 0.0249 


a 


• 

u 


9.6588 

- 0.0278 

0.0000 

- 31.8268 


u 


• 

q 


- 9.9359 

0.0011 

- 8.0298 

0.0836 


q 


• 

e 


0.0000 

0.0000 

1.0000 

0.0000 


e 




- 



— 

i 

- - 






- 0.0360 

- 0.2260 







0.0000 

- 4.5082 


$se 




+ 

- 3.8859 

0.2500 


_ 6f 

J 





0.0000 

0.0000 

1 


y = 

Cx + 

Du 





j 

r ~ 









— 1 

a z 


-215.4679 

-0.3126 - 

16.4373 

3.8543" 


- - 









a 


a 


1.0000 

0.0000 

0.0000 

0.0000 


u 


u 

=r 

0.0000 

1.0000 

0.0000 

0.0000 


q 


q 


0.0000 

0.0000 

1.0000 

0.0000 

1 

0 


e 


0.0000 

0.0000 

0.0000 

1.0000 

















-6.6154 

-41.5528 








0.0000 

0.0000 


1 

1 






0.0000 


5 se 




+ 

0.0000 


fi f 







j 0.0000 

0.0000 








0.0000 

0.0000 

i 
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Table B.12 Cessna 402B Longitudinal Mathematical Model 
For Sea Level Climb (CG = 0.14) 


X 

= Ax 

+ Bu 






• 






-i 


— 


a 


-1.3313 

-0.0014 

0.9163 

-0.0120 


a 


• 

u 


12.7885 

-0.0228 

0.0000 

-32.0688 


u 


• 

q 


-12.2848 

0.0026 

-8.6919 

0.0433 


q 


e 


0.0000 

0.0000 

1.0000 

0.0000 


0 







- 


- - 





— 

- 






-0.0389 

-0.2593 


~ i 




0.0000 

-5.9354 


$se 



+ 









-4.8180 

0.2727 


_$f 





0.0000 

0.0000 

4 



y = 

Cx + Du 




a z 


-280.7059 

-0.2952 ■ 

-17.6482 

2.8614" 

a 


1.0000 

0.0000 

0.0000 

0.0000 

u 


0.0000 

1.0000 

0.0000 

0.0000 

q 


0.0000 

0.0000 

1.0000 

0.0000 

e 


0.0000 

0.0000 

0.0000 

1.0000_ 




"-8.1937 

-54.6737" 




0.0000 

0.0000 



1 

+ 

0.0000 

0.0000 




0.0000 

0.0000 




_ 0.0000 

0.0000_ 



5 se 

6 f 
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Table B.13 Cessna 402B Longitudinal Mathematical Model 
For Climb at 5000 Feet (CG = 0.14) 


x = Ax + Bu 


• 


a 


« 


u 


• 

q 

- 

• 


0 


“ 



-1.2403 

12.5650 

■13.0067 

0.0000 


- 0.0012 

- 0.0212 

0.0019 

0.0000 


0.9282 

0.0000 

-8.0812 

1.0000 

■0.0361 

0.0000 

■4.8260 

0.0000 


-0.0105 

-32.0808 

0.0352| 

0.0000 

-0.2415 

-5.9355 

0.1403 

0.0000 


a 
u 

q 

e 

, J 

$se 


y = Cx + Du 


a z 


“-281.9648 

-0.2728 

-16.3227 

2.8603 



a 


1.0000 

0.0000 

0.0000 

0.0000 


vX 

u 


0.0000 

1.0000 

0.0000 

0.0000 


u 

q 


0.0000 

0.0000 

1.0000 

0.0000 

| 

q 

e 

L — 

e 


0.0000 

0.0000 

0.0000 

1.0000_ 



•8.2114 

0.0000 

0.0000 

0.0000 

0.0000 


-54.9016 

0.0000 

0.0000 

0.0000 

0.0000 


5 se| 
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Table B.14 Cessna 402B Longitudinal Mathematical Model 
For Cruise at 20,000 Feet (CG = 0.14) 


x = Ax + Bu 


-1.2336 

18.3366 

-20.9806 

0.0000 


■0.0005 

■0.0178 

0.0009 

0.0000 


0.9569 

0.0000 

-7.6766 

1.0000 

■0.0342 

0.0000 

-7.1985 

0.0000 


0.0000 

-32.1740 

0.0000 

0.0000 

-0.2422 

-9.2750 

0.2869 

0.0000 


“-441.5215 

-0.1790 

-15.4261 

1.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 


15.4261 1.5492 

0.0000 0.0000 

0.0000 0.0000 

1.0000 0.0000 

0.0000 1.0000 

-12.2549 -86.6865 

0.0000 0.0000 

0.0000 0.0000 

0.0000 0.0000 

0.0000 0.0000 
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Table B.15 Cessna 402B Longitudinal Mathematical Model 
For Sea Level Approach (CG » 0.14) 


x = Ax + Bu 


• 

a 


- 1.0446 

- 0.0024 

0.9166 

0.0102 


— . — , 

a 

• 

u 


21.0420 

- 0.0525 

0.0000 

- 32.1299 


u 

• 

q 


- 7.1555 

0.0053 

- 6.6358 

- 0.0281 


q 

• 

e 


0.0000 

0.0000 

1.0000 

0.0000 


e 



- 



— 


. - 





— 

- 






- 0.0297 

- 0.0569 






0.0000 

- 3.1464 


^se 




+ 








- 2.7946 

0.8094 







0.0000 

0.0000 


y = 

Cx + 

Du 






a z 


-167.3940 

-0.3846 - 

-13.3646 

2.7529 


— n 









a 


a 


1.0000 

0.0000 

0.0000 

0.0000 


u 


u 

= 

0.0000 

1.0000 

0.0000 

0.0000 


q 


q 


0.0000 

0.0000 

1.0000 

0.0000 


e 


e 


0.0000 

0.0000 

0.0000 

1.0000_ 


^ 


— — 












"-4 .7529 

-9.118l“ 








0.0000 

0.0000 


r 









$se 




+ 

0.0000 

0.0000 








0.0000 

0.0000 


— 

— 





0.0000 

0 . 0000 _ 
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Table B.16 Cessna 402B Lateral-Directional 

Mathematical Model For Sea Level Take-off 


x = Ax + Bu 


• 

B 


-0.1547 

0.1178 

• 

P 


-2.9322 

-2.4155 

• 

r 


2.5862 

-0.3308 

• 

$ 


0.0000 

1.0000 


-0.9939 0.1688 


"•1 

0.3692 -0.0062 


p ! 

-0.3206 -0.0063 


r i 

t 

0.2716 0.0000 


<t> ; 

— 


. - 

0.0000 0.0139 


-1.9624 0.2521 


$df 

-0.0615 -0.5288 


«sr 

0.0000 0.0000 



y = Cx + Du 


a y 


~-28 . 4455 

21.6543 

1.1186 

-0.0173’ 

B 


1.0000 

0.0000 

0.0000 

0.0000 

p 

= 

0.0000 

1.0000 

0.0000 

0.0000 

r 


0.0000 

0.0000 

1.0000 

0.0000 

4> 


0.0000 

L 

0.0000 

0.0000 

1.0000 



' 

*0.0000 

2.5606 




0.0000 

0.0000 



+ 

0.0000 

0.0000 




0.0000 

0.0000 




0.0000 

o.oooo_ 


B 

P 

r 

<t> 


5 df 

5 sr 
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Table B.17 Cessna 402B Lateral-Directional 

Mathematical Model For Sea Level Climb 


x = Ax + Bu 


— "1 

• 






8 


-0.1879 

0.0874 

-0.9971 

0.1505 

• 

P 


-3.7107 

-2.6275 

0.3918 

-0.0070 

• 

r 


3.7138 

-0.2901 

-0.3503 

-0.0065 

• 


0.0000 

1.0000 

0.1700 

0.0000 



_ 



- 


0.0000 

-2.6247 

-0.0611 

0.0000 


0.0162 

0.3362 

-0.7013 

0.0000 


6 df 

«sr 


y = Cx + Du 


a y 


-39.6290 

18.4390 

0.6125 

0.0193 

8 


1.0000 

0.0000 

0.0000 

0.0000 

p 

= 

0.0000 

1.0000 

0.0000 

0.0000 

r 


0.0000 

0.0000 

1.0000 

0.0000 

4>_ 


0.0000 

0.0000 

0.0000 

1.0000_ 



] 

“o.oooo 

3.4133 




0.0000 

0.0000 



+ 

0.0000 

0.0000 




0.0000 

0.0000 




0.0000 

o.oooo_ 


Table B.18 Cessna 402B Lateral-Directional 

Mathematical Model For Climb at 5000 Feet 


x = Ax + Bu 


- ' 


— 





— - 


B 


-0.1742 

0.0876 

-0.9969 

0.1397 


B 


• 

P 


-3.7132 

-2.4327 

0.3639 

-0.0060 


P 


• 

r 


3.7107 

-0.2692 

-0.3246 

-0.0056 


r 


i 


0.0000 

1.0000 

0.1653 

0.0000 


<t> 



i 

- 



— 


- — 





0.0000 

0.0150 





-2.6244 

0.3363 


$df 



+ 








-0.0612 

-0.7012 


$sr 




0.0000 

0.0000 



y = Cx + Du 


a y" 


“-39.6047 

19.9233 

0.6936 

0.0113 

B 


1.0000 

0.0000 

0.0000 

0.0000 

P 

= 

0.0000 

1.0000 

0.0000 

0.0000 

r 


0.0000 

0.0000 

1.0000 

0.0000 

_ 4 > 


0.0000 

q.oooo 

0.0000 

1.0000_ 



1 

“o.oooo 

3 . 4121 " 




0.0000 

0.0000 



+ 

0.0000 

0.0000 




0.0000 

0.0000 




0.0000 

0.0000 



$df 

5 sr 
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Table B.19 Cessna 402B Lateral-Directional 

Mathematical Model For Cruise at 20,000 Feet 


x = Ax + Bu 


8 

P 

r 

$ 


-0.1843 

-5.3309 

6.3295 

0.0000 


0.0482 

-2.3284 

■0.1677 

1.0000 


y = Cx + Du 


1 

>1 

fO 


8 


P 

= 

r 


_4> _ 

- 


-0.9993 

0.3048 

-0.3141 

0.0482 

0.0000 

-4.1555 

-0.0435 

0.0000 


0.0898 

-0.0039 

■0.0032 

0.0000 

0.0151 

0.5331 

-1.1090 

0.0000 


- 

8 


P 


r 

-J 

L 


$df 

®sr 


-65.9673 

17.2417 

1.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

q.oooo 


0.2434 

0.0000 

0.0000 

1.0000 

0.0000 

" 0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


0.0107 

0.0000 

0.0000 

0.0000 

1 . 0000 _ 

5.4192 

0.0000 

0.0000 

0.0000 

0.0000 


8 

p 

r 

♦ 


$df 


»sr 
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Table B.20 Cessna 402B Lateral-Directional 

Mathematical Model For Sea Level Approach 


-0.1452 

-2.1765 

2.1817 

0.0000 


0.0871 

■2.0130 

-0.2224 

1.0000 


■0.9971 

0.3034 

-0.2692 

0.0348 

0.0000 

-1.5409 

-0.0357 

0.0000 


0.2008 

■0.0072 

-0.0067 

0.0000 

0.0125 

0.1974 

-0.4117 

0.0000 


-23.2752 

13.9586 

0.4583 

1.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

q.oooo 

0.0000 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


0.0216 

0.0000 

0.0000 

0.0000 

1 . 0000 _ 

2.0036 

0.0000 

0.0000 

0.0000 

0.0000 
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Max . 

RMS 

Max. 

a z 

(f/s2) 

16.2 

5.7 

10.3 

a 

(deg) 

4.4 

1.6 

5.4 

u 

(f/s) 

2.1 

1.2 

3.4 

q 

(deg/s) 

2.4 

1.0 

3.7 

e 

(deg) 

2.6 

0.9 

3.5 

fi se 

(deg) 

- 

- 

4.4 

$f 

(deg) 

- 

- 

15.0 

^se 

(deg/s) 

- 

- 

7.9 

$f 

(deg/s) 

- 

- 

83.8 

* 


B.23 


inal RQAS Response 

e-off — Control-Rate-Weighting 


seline 

RQAS 


Percent 
Reduct ion 


RMS 

Max . 

RMS 

3.4 

36% 

40% 

1.8 

-23% 

-12% 

1.9 

-62% 

-58% 

1.6 

-54% 

-60% 

1.3 

-35% 

-44% 

1.6 

_ 

_ 


5.6 


Table B.21 Summary of Longitudinal RQAS Response 

for a Sea-Level take-off — Control-Rate-Weighting 
(Continued) 


System Eigenvalues : 


Open Loop 


Baseline RQAS 


Frequency Damping 


Frequency Damping 


Phugoid 

0.15 

0.02 

0.09 

0.87 

Phugoid 

0.15 

0.02 

0.09 

0.87 

Short Period 

2.1 

1.0 

6.3 

0.70 

Short Period 

6.6 

1.0 

6.3 

0.70 

Servo 

10.0 

1.0 

2.4 

1.0 

Servo 

10.0 

1.0 

11.6 

1.0 

Filter 

- 

- 

14.6 

1.0 

Filter 

- 

- 

36.6 

1.0 


Note 1: Baseline RQAS Frequency and Damping are based on 

roots of the W' -Plane. 


Note 2: Control-Rate-Weighting introduces 2 filter states. 


Optimal Gain Matrix : 


K 


5.4211 -0.0021 -4.0595 

110.92 0.1610 12.198 


-9.0508 

22.167 


7.9180 

1.5731 


1.9853 

26.910 


Weighting Matrices — Diagonal Elements : 


Q’ = (0.05, 10.0, 
R' = (0.77, 3.06) 
S’ = (0.07, 0.06) 


0.0001, 15.0, 20.0) 
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Table B. 22 Summary of Longitudinal RQAS Response 

for a Sea-Level Climb — Control-Rate-Weighting 


Time Domain Response: 


Open 

Loop 


Baseline 

RQAS 


Percent 

Reduction 



Max. 

RMS 

Max. 

RMS 

Max. 

RMS 

a z 

( f/s2) 

13.1 

6.3 

12.2 

3.4 

7% 

46% 

a 

(deg) 

2.8 

1.3 

3.6 

1.5 

-29% 

-15% 

u 

(f/s) 

1.9 

1.1 

4.1 

2.4 

-116% 

-118% 

q 

(deg/s) 

1.8 

0.9 

2.7 

1.4 

-50% 

-56% 

6 

(deg) 

1.6 

1.0 

2.0 

1.1 

-25% 

-10% 

$se 

(deg) 

- 

- 

4.3 

2.3 

- 

- 

$f 

(deg) 

- 

- 

11.8 

3.4 

- 

- 

®se 

(deg/s) 

- 

- 

11.4 

- 

- 

- 

$f 

(deg/s) 

- 

- 

91.1 

- 

- 

- 
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Table B.22 Summary of Longitudinal RQAS Response 

for a Sea-Level Climb — Control-Rate-Weighting 
(Continued) 


System Eigenvalues : 


Open Loop 
Frequency Damping 


Baseline RQAS 
Frequency Damping 


Phugoid 

Phugoid 

0.15 

0.04 

0.09 

0.79 

0.15 

0.04 

0.09 

0.79 

Short Period 

2.4 

1.0 

6.0 

0 . 68 

Short Period 

7.1 

1.0 

6.0 

0.68 

Servo 

10.0 

1.0 

3.5 

1.0 

Servo 

10.0 

1.0 

11.8 

1.0 

Filter 

- 

- 

15.3 

1.0 

Filter 

- 

— 

63.8 

1.0 

Note 1: Baseline 

roots of 

! RQAS Frequency 
the W' -Plane. 

and Damping are 

based on 

Note 2: Control- 

Rate- 

■Weighting 

introduces 2 filter state: 


Matrix: 


Optimal 

K = 


Gain 
27.930 
171.39 


0.0123 

0.1719 


-2.5093 

14.655 


-13.634 

24.772 


8.8770 

3.0441 


6.1749 

38.537 


Weighting Matrices 

Q' = (0.05, 10.0, 
R- = (0.12, 0.33) 
S' = (0.02, 0.03) 


— Diagonal Elements : 
0.0001, 3.0, 20.0) 
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Table B.23 Summary of Longitudinal RQAS Response 

for a Climb at 5000 feet— Control-Rate-Weighting 


Time Domain Response : 

Open 

Loop 


Baseline 

RQAS 


Percent 

Reduction 



Max. 

RMS 

Max. 

RMS 

Max. 

RMS 

a z 

(f/s2) 

10.2 

3.6 

7.3 

2.1 

28% 

42% 

a 

(deg) 

2.2 

0.8 

3.1 

1.2 

-41% 

-50% 

u 

(f/s) 

2.4 

1.4 

4.2 

2.4 

-75% 

-71% 

q 

(deg/s) 

1.5 

0.7 

.2.1 

0.9 

-40% 

-28% 

6 

(deg) 

2.2 

1.1 

2.4 

1.2 

-9% 

-9% 

fi se 

(deg) 

- 

- 

4.2 

2.2 

- 

- 

$f 

(deg) 

- 

- 

9.0 

4.1 

- 

- 

$se 

(deg/s) 

- 

- 

6.4 

- 

- 

- 

«f 

(deg/s) 

- 

- 

35.3 

- 

- 

— 
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Table B.23 Summary of Longitudinal RQAS Response 

for a Climb at 5000 feet — Control-Rate-Weighting 
(Continued) 


System Eigenvalues: 


Open Loop Baseline RQAS 


Frequency Damping Frequency Damping 


Phugoid 

Phugoid 

0.14 

0.04 

0.08 

0.80 

0.14 

0.04 

0.08 

0.80 

Short Period 

2.5 

1.0 

6 . 1 

0.63 

Short Period 

6.2 

1.0 

6.1 

0.63 

Servo 

10.0 

1.0 

3.5 

1.0 

Servo 

10.0 

1.0 

11.6 

1.0 

Filter 

- 

- 

15.6 

1.0 

Filter 

- 

- 

64.1 

1.0 


Note 1: Baseline RQAS Frequency and Damping are based on 

roots of the W' -Plane. 

Note 2: Control-Rate-Weighting introduces 2 filter states. 


Matrix: 


Optimal 
K = 


Gain 
25.662 
171.79 


0.0105 

0.1520 


-3.2351 

14.182 


-12.532 

25.136 


9.2575 

3.0844 


6.1700 

38.631 


Weighting Matrices — Diagonal Elements : 


Q’ 

« (0.05, 

15.0, 

R' 

* (0.12, 

0.33) 

S’ 

= (0.02, 

0.03) 


0.0001, 3.0, 20.0) 
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Max. 

RMS 

Max. 

a Z 

(f/s2) 

7.0 

2.9 

4.0 

a 

(deg) 

0.9 

0.4 

1.4 

u 

(f/s) 

1.5 

0.9 

2.4 

q 

(deg/s) 

1.2 

0.6 

1.3 

0 

(deg) 

1.3 

0.7 

1.4 

«se 

(deg) 

- 

- 

2.0 

5 f 

(deg) 

- 

- 

4.9 

6 g0 

(deg/s) 

- 

- 

7.0 

Sf 

(deg/s) 

- 

- 

25.3 
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nal RQAS Response 
0 feet — Control-Rate-Weighting 


ieline 

IQAS 


Percent 

Reduction 


RMS 

Max. 

RMS 

1.5 

43% 

48% 

0.6 

-56% 

-50% 

1.5 

-60% 

-67% 

0.6 

-8% 

0% 

00 

• 

o 

-8% 

-14% 

0.9 

- 

- 

2.3 

- 

- 





Table B.24 Summary of Longitudinal RQAS Response 

for a Cruise at 20000 ft — Control-Rate-Weighting 
(Continued) 


System Eigenvalues ; 



Open 

Loop 

Baseline 

RQAS 


Freguency Damping 

Frequency 

Damping 

Phugoid 

0.10 

0.08 

0.08 

0.78 

Phugoid 

0.10 

0.08 

0.08 

0.78 

Short Period 4.5 

0.9 

10.2 

0.48 

Short Period 4.5 

0.9 

10.2 

0.48 

Servo 

10.0 

1.0 

1.3 

1.0 

Servo 

10.0 

1.0 

11.5 

1.0 

Filter 

- 

- 

27.9 

1.0 

Filter 

— 

— 

125 

1.0 

Note 1: 

Baseline RQAS Frequency and 
roots of the W' -Plane. 

Damping are based on 


Note 2: Control-Rate-Weighting introduces 2 filter states. 


Optimal 

Gain Matrix: 





K = 

29.784 

-0.0443 

-15.303 

-5.6702 

18.406 

4.8751 


254.71 

0.0959 

15.162 

21.686 

3.9323 

55.227 


— 





M, 


Weighting Matrices — Diagonal Elements : 


Q' 

= (0.05, 

10.0, 0.0001, 15.0, 20.0) 

R' 

= (0.06, 

0.11) 

S' 

= (0.01, 

0.01) 
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Table B.25 Summary of Longitudinal RQAS Response 
for an Approach at Sea-Level 
Control-Rate-Weighting 


Time Domain Response : 

Open Baseline Percent 

Loop RQAS Reduction 



Max. 

RMS 

Max. 

RMS 

Max. 

RMS 

a z 

(f/s2) 

12.9 

4.9 

10.7 

3.9 

17% 

20% 

a 

(deg) 

4.8 

1.8 

4.7 

1.7 

2% 

6% 

u 

(i/s) 

2.9 

1.6 

3.5 

2.1 

-21% 

-31% 

q 

(deg/s) 

2.2 

1.0 

3.5 

1.7 

-59% 

-70% 

e 

(deg) 

2.7 

1.3 

3.5 

1.5 

-30% 

-15% 

$se 

(deg) 

- 

- 

4.4 

1.6 

— 


«f 

(deg) 

- 

- 

15.5 

6.4 

— 

“ 

6 gg 

(deg/s) 

- 

- 

16.8 

— 

— 


*f 

(deg/s) 

- 

- 

101.4 

* 

- 

— 

— 
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Table B.25 Summary of Longitudinal RQAS Response 
for an Approach at Sea-Level 
Control-Rat e-Weighting (Continued) 


System Eigenvalues ; 



Open 

Loop 

Baseline 

RQAS 


Frequency Damping 

Frequency 

Damping 

Phugoid 

0.21 

0.15 

0.19 

0.93 

Phugoid 

0.21 

0.15 

0.19 

0.93 

Short Period 1.8 

1.0 

7.4 

0.67 

Short Period 5.4 

1.0 

7.4 

0.67 

Servo 

10.0 

1.0 

3.2 

1.0 

Servo 

10.0 

1.0 

10.6 

1.0 

Filter 

- 

- 

21.4 

1.0 

Filter 

— 

- 

32.4 

1.0 

Note Is 

Baseline RQAS Frequency 
roots of the W' -Plane. 

and Damping are based on 


Note 2: Control-Rate-Weighting introduces 2 filter states. 


Optimal Gain Matrix : 


K 


-12.514 

138.74 


-0.7461 

0.3300 


-14.385 

31.573 


-23.176 

74.169 


16.034 

-1.3032 


-0.1902 

23.9256 


Weighting Matrices — Diagonal Elements : 


Q' 

= (0.01, 

10.0, 

R' 

= (2.01, 

1.002 

S' 

= (0.01, 

0.002 


0.0001, 15.0, 20.0) 
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Open Baseline Percent 

Loop RQAS Reduction 



Max. 

RMS 

Max. 

RMS 

Max. 

RMS 

a y 

( f/s2 ) 

5.2 

2.2 

3.4 

1.1 

35% 

50% 

B 

(deg) 

7.6 

2.9 

5.5 

2.2 

28% 

24% 

P 

(deg/s) 

5.6 

2.5 

1.7 

0.5 

70% 

8% 

r 

(deg/s) 

11.3 

4.5 

6.0 

2.7 

47% 

40% 

♦ 

(deg) 

3.7 

1.6 

1.8 

0.8 

51% 

50% 

«df 

(deg) 

- 

- 

8.5 

2.7 

- 

- 

«sr 

(deg) 

- 

- 

13.4 

6.0 

- 

- 

, $df % 
(deg/s) 

- 

- 

43.4 

- 

- 

- 

$sr 

(deg/s) 

- 

- 

32.8 

- 

- 

- 
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Table B.26 Summary of Lateral-Directional RQAS Response 

for a Take-off at Sea-Level — Control-Rate-Weighting 
(Continued) 


System Eigenvalues : 

Open Loop Baseline RQAS 

Frequency Damping Frequency Damping 


Spiral Mode 

0.04 

-1.0 

0.36 

1.0 

Roll Mode 

2.5 

1.0 

2.8 

1.0 

Dutch Roll 

1.8 

0.13 

2.0 

0.55 

Dutch Roll 

1.8 

0.13 

2.0 

0.55 

Servo 

10.0 

1.0 

5.9 

0.85 

Servo 

10.0 

1.0 

5.9 

0.85 

Filter 

- 

- 

12.6 

1.0 

Filter 

- 

— 

21.3 

1.0 

Note 1: Baseline RQAS 

roots of the W 

Frequency 
'' -Plane. 

and Damping are 

based on 


Note 2: Control-Rate-Weighting introduces 2 filter states. 


Optimal Gain Matrix: 


K 


29.886 

-4.5862 


-17.330 

3.1201 


-7.8019 

-9.9103 


-13.481 

2.4044 


19.073 

-0.2559 


-1.5378 

5.3043 


Weighting Matrices — Diagonal Elements : 


Q' = (0.025, 

o 

o 

• 

r- 

R' = (3.007, 

1.07) 

S' = (0.007, 

0.07) 


0.07, 


5.0, 1.5) 
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Open Baseline Percent 

Loop RQAS Reduction 



Max. 

RMS 

Max. 

RMS 

Max. 

RMS 

a y 

( f/s2 ) 

6.1 

2.7 

3.5 

1.3 

43% 

52% 

B 

(deg) 

6.2 

2.8 

5.8 

1.9 

6% 

32% 

P 

(deg/s) 

6.4 

2.7 

1.5 

0.5 

77% 

81% 

r 

(deg/s) 

10.6 

4.8 

7.1 

2.8 

33% 

42% 

^ ♦ 

(deg) 

3.2 

1.3 

1.0 

0.4 

69% 

69% 

«df 

1 (deg) 

- 

- 

4.8 

2.1 

- 

— 

«sr 

(deg) 

- 

- 

14.7 

5.8 

— 


$df 

(deg/s) 

- 

- 

30.0 

- 

— 

— 

$sr 

(deg/s) 

- 

- 

30.8 

- 

- 

— 
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Table B.27 Summary of Lateral-Directional RQAS Response 

for a Sea-Level Climb — Control-Rate-Weighting 
(Continued) 


System Eigenvalues : 


Open Loop Baseline RQAS 


Frequency Damping Frequency Damping 


Spiral Mode 

0.02 

-1.0 

0.59 

1.0 

Roll Mode 

2.7 

1.0 

3.3 

1.0 

Dutch Roll 

2.1 

0.12 

2.8 

0.51 

Dutch Roll 

2.1 

0.12 

2.8 

0.51 

Servo 

10.0 

1.0 

5.5 

0.88 

Servo 

10.0 

1.0 

5.5 

0.88 

Filter 

- 

- 

12.6 

1.0 

Filter 

- 

— 

21.3 

1.0 

Note 1: Baseline 

roots of 

RQAS Frequency 
the W' -Plane. 

and Damping are based on 

Note 2: Control- 

Rate- 

-Weighting 

introduces 2 filter 

states 


Optimal Gain Matrix: 


K 


22.429 

-1.9138 


-9.4558 

3.7262 


-6.6832 

-14.873 


-12.139 

3.5634 


16.148 

-0.4145 


-0.6806 

7.2291 


Weighting Matrices — Diagonal Elements : 

Q' = (0.01, 15.0, 0.07, 0.5, 1.5) 

R' * (2.007, 0.52) 

S' = (0.007, 0.02) 
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irectional RQAS Response 
ft — Control-Rate-Weighting 


seline 

*QAS 


Percent 

Reduction 


RMS 

Max. 

RMS 

0.9 

28% 

25% 

1.2 

3% 

8% 

0.4 

65% 

64% 

1.6 

-3% 

16% 

0.2 

55% 

60% 

1.6 

- 

- 

3.2 


_ 


Table B.28 Summary of Lateral-Directional RQAS Response 

for a Climb at 5000 ft — Control-Rate-Weighting 
(Continued) 


System Eigenvalues : 



Open 

Loop 

Baseline 

RQAS 


Freguency Damping 

Frequency 

Dampint 

Spiral 

Mode 0.02 

-1.0 

0.58 

1.0 

Roll Mode 2.5 

1.0 

3.2 

1.0 

Dutch ] 

Roll 2.1 

0.11 

2.7 

0.50 

Dutch ] 

Roll 2.1 

0.11 

2.7 

0.50 

Servo 

10.0 

1.0 

5.6 

0.85 

Servo 

10.0 

1.0 

5.6 

0.85 

Filter 

- 

- 

12.6 

1.0 

Filter 

- 

— 

21.8 

1.0 

Note 1 

: Baseline RQAS 

Frequency 

and Damping are based on 


roots of the 

W' -Plane. 



Note 2 

: Control-Rate- 

Weighting 

introduces 2 filter 

states 


Opt 

K 


imal Gain Matrix : 

25.050 -10.715 

-2.1522 3.9191 


-7.7077 

-15.067 


-11.876 

3.7219 


16.368 

-0.4313 


-0.6975 

7.2545 


Weighting Matrices — Diagonal Elements : 

Q' = (0.01, 15.0, 0.07, 0.5, 1.5) 

R’ = (2.007, 0.52) 

S' * (0.007, 0.02) 
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Table B.29 Summary of Lateral-Directional RQAS Response 

for a Cruise at 20000 ft — Control-Rate-Weighting 


Time Domain Response : 

Open Baseline Percent 

Loop RQAS Reduction 



Max. 

RMS 

Max. 

RMS 

Max. 

RMS 

a y 

( f/s2 ) 

2.4 

1.1 

1.7 

0.5 

29% 

55% 

S 

(deg) 

1.9 

0.8 

1.5 

0.5 

21% 

38% 

P 

(deg/s) 

2.3 

1.0 

0.8 

0.3 

65% 

70% 

r 

(deg/s) 

4.1 

1.8 

2.3 

0.9 

44% 

50% 

4> 

(deg) 

1.0 

0.4 

0.4 

0.1 

60% 

25% 

, fi df 

(deg) 

- 

- 

2.0 

0.8 

- 

- 

, 5 sr 

(deg) 

- 

- 

4.5 

1.8 

- 

- 

, ® df X 

(deg/s) 

- 

- 

13.8 

- 

- 

- 

$sr 

(deg/s) 

- 

- 

13.4 

- 

- 

- 
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Table B.29 Summary of Lateral-Directional RQAS Response 

for a Cruise at 20000 ft — Control-Rate-Weighting 
(Continued) 


System Eigenvalues : 



Open 

Loop 

Baseline 

RQAS 


Frequency Damping 

Frequency 

Damping 

Spiral Mode 

0.003 

-1.0 

0.80 

1.0 

Roll Mode 

2.4 

1.0 

2.6 

1.0 

Dutch Roll 

2.6 

0.08 

3.5 

0.45 

Dutch Roll 

2.6 

0.08 

3.5 

0.45 

Servo 

10.0 

1.0 

6.3 

0.72 

Servo 

10.0 

1.0 

6.3 

0.72 

Filter 

- 

- 

14.1 

1.0 

Filter 

- 

- 

23.6 

1.0 

Note 1: Baseline RQAS 

Frequency 
W' -Plane. 

and Damping are based on 

roots 

of the 



Note 2: Control-Rate- 

Weighting 

introduces 2 filter 

states. 


Matrix: 


Optimal 
K = 


Gain 
41.555 
11.781 


-10.172 

3.6869 


-11.371 

-14.791 


-10.074 

4.9004 


17.282 

-0.5793 


-0.7426 

8.5732 


Weighting Matrices — Diagonal Elements : 

Q’ = (0.01 ,15.00, 0.07, 0.5, 1.5) 

R’ * (2.007, 0.52) 

S' = (0.007, 0.02) 
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Table B.30 Summary of Lateral- 
for an Approach at 



Open 

Loop 



Max. 

RMS 

Max. 

a y 

( f /s2 ) 

5.8 

2.4 

2.6 

0 

(deg) 

10.2 

4.2 

5.6 

P 

(deg/s) 

7.5 

3.1 

2.0 

r 

(deg/s) 

10.5 

5.8 

6.4 

* 

(deg) 

3.8 

2.0 

1.1 

, 5 df % 
(deg) 

- 

- 

4.6 

$sr 

(deg) 

- 

- 

11.7 

«df 

(deg/s) 

- 

- 

14.3 

5sr 

(deg/s) 

- 

- 

21.6 
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rectional RQAS Response 
a-Level — Control-Rate-Weighting 


seline 

RQAS 


Percent 

Reduction 


RMS 

Max. 

RMS 

1.2 

55% 

50% 

2.5 

45% 

40% 

0.9 

73% 

71% 

3.4 

39% 

41% 

0.6 

71% 

70% 

2.1 

- 

- 

6.2 

- 

- 





Table B.30 Summary of Lateral-Directional RQAS Response 

for an Approach at Sea-Level — Control-Rate-Weighting 
(Continued) 


System Eigenvalues : 

Open Loop 
Frequency Damping 


Spiral Mode 

0.003 

-1.0 

Roll Mode 

2.1 

1.0 

Dutch Roll 

1.6 

0.11 

Dutch Roll 

1.6 

0.11 

Servo 

10.0 

1.0 

Servo 

10.0 

1.0 

Filter 

- 

- 

Filter 

- 

- 


Baseline RQAS 
Frequency Damping 


0.55 

1.0 

2.7 

1.0 

1.8 

0.40 

1.8 

0.40 

5.1 

1.00 

7.6 

1.00 

11.9 

1.0 

19.4 

1.0 


Note 1: Baseline RQAS Frequency and Damping are based on 

roots of the W' -Plane. 


Note 2: Control-Rate-Weighting introduces 2 filter states. 


Optimal Gain Matrix : 

10.580 -8.5909 

1.9229 2.8267 


K ■ 


-1.1162 

-14.721 


-12.695 

2.9242 


15.332 

-0.2052 


-0.4045 

7.7058 


Weighting Matrices — Diagonal Elements : 

Q' = (0.01 ,15.00, 0.07, 0.5, 1.5) 

R* - (2.007, 1.02) 

S' = (0.007, 0.02) 
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Max. 

RMS 

Max. 

a z 

( f/s2) 

16.2 

5.7 

10.4 

a 

(deg) 

4.4 

1.6 

5.4 

u 

(f/s) 

2.1 

1.2 

3.1 

q 

(deg/s) 

2.4 

1.0 

3.4 

6 

(deg) 

2.6 

0.9 

3.3 

5 se 

(deg) 

- 

- 

4.9 

, 6f 

(deg) 

- 

- 

15.0 

6 se 

(deg/s) 

- 

- 

9.4 

(deg/s) 

- 

- 

103.9 
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inal RQAS Response 
Level — Output-Weighting 


seline 

RQAS 


Percent 

Reduction 


RMS 

Max. 

RMS 

3.3 

36% 

42% 

CO 

• 

r-f 

-23% 

-13% 

1.8 

-48% 

-50% 

1.5 

-42% 

-50% 

1.2 

-27% 

-33% 

1.8 

- 

- 

5.6 





Table B.31 Summary of Longitudinal RQAS Response 

for Take-Off at Sea-Level — Output-Weighting 
(Continued) 


System Eigenvalues : 

Open Loop 
Frequency Damping 


Baseline RQAS 
Frequency Damping 


Phugoid 

Phugoid 

0.15 

0.02 

0.09 

1.00 

0.15 

0.02 

0.11 

1.00 

Short Period 

2.1 

1.0 

8.7 

0.87 

Short Period 

6.6 

1.0 

8.7 

0.87 

Servo 

10.0 

1.0 

2.1 

1.0 

Servo 

10.0 

1.0 

11.1 

1.0 


Note 1: Baseline RQAS Frequency and Damping are based on 

roots of the W' -Plane. 


Opt imal 
K = 


Gain Matrix : 
-0.1664 -0.0011 

4.1356 0.0060 


-0.6194 

0.4927 


-1.3505 

1.1368 


Weighting Matrices — Diagonal Elements : 


Q' = (0.04, 1.00, 0.0001, 18.0, 40.0) 
R' = (8.0, 2.5) 
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Table B.32 Summary of Longitudinal RQAS Response 

for a Sea-Level Climb — Output-Weighting 


Time Domain Response : 

Open 

Loop 


Baseline 

RQAS 


Percent 

Reduction 



Max. 

RMS 

Max. 

RMS 

Max. 

RMS 

a z 

(f/s2) 

13.1 

6.3 

11.8 

3.3 

10% 

48% 

a 

(deg) 

2.8 

1.3 

3.7 

1.6 

-32% 

-23% 

u 

(f/s) 

1.9 

1.1 

4.0 

2.3 

-111% 

-110% 

q 

(deg/s) 

1.8 

0.9 

2.5 

1.2 

-39% 

-33% 

0 

(deg) 

1.6 

1.0 

1.8 

1.1 

-13% 

-10% 

^se 

(deg) 

- 

- 

5.0 

2.5 

- 

- 

«f 

(deg) 

- 

- 

11.7 

5.4 

- 

- 

$se 

(deg/s) 

- 

- 

22.3 

- 

- 

- 

«f 

(deg/s) 

- 

- 

98.8 

- 

- 

- 
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Table B.32 Summary of Longitudinal RQAS Response 

for a Sea-Level Climb — Output-Weighting 
(Continued) 


System Eigenvalues : 


Open Loop 


Baseline RQAS 


Frequency Damping 


Frequency Damping 


Phugoid 

0.15 

0.04 

0.09 

0.89 

Phugoid 

0.15 

0.04 

0.09 

0.89 

Short Period 

2.4 

1.0 

11.6 

0.66 

Short Period 

7.1 

1.0 

11.6 

0.66 

Servo 

10.0 

1.0 

1.9 

1.0 

Servo 

10.0 

1.0 

11.3 

1.0 


Note 1: Baseline RQAS Frequency and Damping are based on 

roots of the W' -Plane. 


Optimal Gain Matrix: 



— 



— 


-0.5014 

-0.0018 

-1.5609 

-2.5312 

K = 

4.5062 

0.0046 

0.5469 

0.9847 


Weighting Matrices — Diagonal Elements : 


Q* = (0.05, 10.0, 0.0001, 18.0, 40.0) 
R’ = (2.2, 1.0) 
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Table B.33 Summary of Longitudinal 
for a Climb at 5000 ft 


Time Domain Response: 


RQAS Response 
- Output-Weighting 


Open Baseline Percent 

Loop RQAS Reduction 



Max. 

RMS 

Max. 

RMS 

Max . 

RMS 

a z 

(f/s2) 

10.2 

3.6 

7.5 

2.1 

26% 

42% 

a 

(deg) 

2.2 

0.8 

3.3 

1.3 

-50% 

-63% 

u 

(f/s) 

2.4 

1.4 

4.1 

2.3 

-71% 

-64% 

q 

(deg/s) 

1.5 

0.7 

1.8 

0.8 

-20% 

-14% 

e 

(deg) 

2.2 

1.1 

2.2 

1.2 

0% 

9% 

«se 

(deg) 

- 

- 

5.0 

2.5 

- 

- 

5 f 

(deg) 

- 

- 

9.6 

2.1 

- 

- 

5se 

(deg/s) 

- 

- 

11.4 

- 

- 

- 

«f 

(deg/s) 

- 

- 

39.5 

- 

- 

- 
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Table B.33 Summary of Longitudinal RQAS Response 

for a Climb at 5000 ft — Output-Weighting 
(Continued) 


System Eigenvalues : 

Open Loop Baseline RQAS 

Frequency Damping Frequency Damping 


Phugoid 

0.14 

0.04 

0.09 

0.89 

Phugoid 

0.14 

0.04 

0.09 

0.89 

Short Period 

2.6 

1.0 

11.5 

0 . 64 

Short Period 

6.2 

1.0 

11.5 

0.64 

Servo 

10.0 

1.0 

2.0 

1.0 

Servo 

10.0 

1.0 

11.2 

1.0 


Note 1: Baseline RQAS Frequency and Damping are based on 

roots of the W' -Plane. 


Optimal Gain Matrix : 


K 


-0.6918 

4.5239 


-0.0016 

0.0041 


-1.6333 

0.5369 


-2.4691 

0.9860 


Weighting Matrices — Diagonal Elements ; 


Q’ = (0.05, 15.0, 0.0001, 18.0, 40.0) 
R’ = (2.2, 1.0) 



B . 48 


Table B.34 Summary of Longitudinal RQAS Response 

for Cruise at 20000 ft — Output-Weighting 


Time Domain Response ; 

Open 

Loop 


Baseline 

RQAS 


Percent 

Reduction 



Max . 

RMS 

Max. 

RMS 

Max. 

RMS 

a z 

(f/s2) 

7.0 

2.9 

4.3 

1.5 

39% 

48% 

a 

(deg) 

0.9 

0.4 

1.4 

0.5 

-56% 

-25% 

u 

(i/s) 

1.5 

0.9 

2.4 

1.5 

-60% 

-67% 

q 

(deg/s) 

1.2 

0.6 

1.3 

0.6 

- 8% 

0% 

6 

(deg) 

1.3 

0.7 

1.4 

0.7 

-8% 

0% 

5 se 

(deg) 

- 

- 

2.3 

0.9 

- 

- 

*f 

(deg) 

- 

- 

4.8 

2.0 

- 

- 

$se 

(deg/s) 

- 

- 

9.0 

- 

- 

- 

$f 

(deg/s) 

- 

- 

23.8 

- 

- 

- 
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Table B.34 Summary of Longitudinal RQAS Response 

for Cruise at 20000 ft — Output-Weighting 
(Continued) 


System Eigenvalues : 


Open Loop 


Baseline RQAS 


Frequency Damping Frequency Damping 


Phugoid 

Phugoid 

0.10 

0.08 

0.08 

0.81 

0.10 

0.08 

0.08 

0.81 

Short Period 

4.5 

0.94 

13.5 

0.51 

Short Period 

4.5 

0.94 

13.5 

0.51 

Servo 

10.0 

1.0 

1.9 

1.0 

Servo 

10.0 

1.0 

11.3 

1.0 


Note 1: Baseline RQAS Frequency and Damping are based on 

roots of the W' -Plane. 


Optimal Gain Matrix : 

-1.2008 -0.0031 -1.7778 -1.5653 

4.7486 0.0020 0.4348 0.6318 


K = 


Weighting Matrices — Diagonal Elements ; 


Q' = (0.05, 10.0, O.OdOl, 18.0, 40.0) 
R' = (2.2, 1.0) 
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Table B.35 Summary of Longitudinal RQAS Response 

for an Approach at Sea-Level — Output-Weighting 


Time Domain Response : 

Open Baseline Percent 

Loop RQAS Reduction 


Max. RMS Max. 


a z 

( f/s2 ) 

12.9 

4.8 

10.6 

a 

(deg) 

4.8 

1.8 

4.5 

u 

(f/s) 

2.9 

1.6 

3.5 

q 

(deg/s) 

2.2 

1.0 

3.9 

e 

(deg) 

2.7 

1.3 

3.8 

$se 

(deg) 

- 

- 

4.8 

s f 

(deg) 

- 

- 

14.6 

5 se 

(deg/s) 

- 

- 

23. 

$f 

(deg/s) 

- 

- 

104. 

ft 


RMS Max . RMS 


3.8 

18% 

21% 

1.7 

- 6% 

- 6% 

2.1 

-21% 

-31% 

1.9 

-77% 

-90% 

1.6 

-41% 

-23% 

1.5 

- 

- 

5.8 



m 


4 

5 
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Table B.35 Summary of Longitudinal RQAS Response 

for an Approach at Sea-Level — Output-Weighting 
(Continued) 


System Eigenvalues : 

Open Loop 
Frequency Damping 


Baseline RQAS 
Frequency Damping 


Phugoid 

0.21 

0.15 

0.18 

0.94 

Phugoid 

0.21 

0.15 

0.18 

0.94 

Short Period 

1.8 

1.00 

6.9 

0.86 

Short Period 

5.4 

1.00 

6.9 

0.86 

Servo 

10.0 

1.0 

4.5 

1.0 

Servo 

10.0 

1.0 

10.3 

1.0 


Note Is Baseline RQAS Frequency and Damping are based on 
roots of the W' -Plane. 


Opt 


K 


imal Gain Matrix : 

-1.1832 -0.0050 

5.3214 0.0135 


-0.7078 

1.0758 


-1.3297 

2.7676 


Weighting Matrices — Diagonal Elements : 


Q' = (0.02, 10.0, O.odoi, 18.0, 40.0) 
R' = (9.5, 3.2) 
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Table B. 36 Summary of Lateral-Directional RQAS Response 
for Take-Off at Sea-Level- — Output-Weighting 


Time Domain Response : 

Open 

Loop 


Baseline 

RQAS 


Percent 

Reduction 



Max. 

RMS 

Max. 

RMS 

Max. 

RMS 

a y 

(f/s2) 

5.2 

2.2 

3.1 

1.0 

40% 

55% 

6 

(deg) 

7.6 

2.9 

5.4 

2.2 

29% 

24% 

P 

(deg/s) 

5.6 

2.4 

1.6 

0.5 

71% 

79% 

r 

(deg/s) 

11.3 

4.5 

5.9 

2.7 

48% 

40% 

<l> 

(deg) 

3.7 

1.6 

1.6 

0.8 

57% 

50% 

, 6 df 

(deg) 

- 

- 

9.7 

2.9 

- 

- 

«sr 

(deg) 

- 

- 

14.4 

6.4 

- 

- 

«df v 
(deg/s) 

- 

- 

55.1 

- 

- 

- 

, ^ Sr V 

(deg/s) 

- 

- 

42.2 

- 

- 

- 
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Table B.36 Summary of Lateral-Directional RQAS Response 
for Take-Off at Sea-Level — Output-Weighting 
(Continued) 


System Eigenvalues : 



Open 

Loop 

Baseline 

RQAS 


Frequency Damping 

Frequency 

Damp i : 

Spiral Mode 

0.04 

-1.0 

0.67 

1.0 

Roll Mode 

2.5 

1.0 

7.5 

1.0 

Dutch Roll 

1.8 

0.13 

2.0 

0.60 

Dutch Roll 

1.8 

0.13 

2.0 

0 . 60 

Servo 

10.0 

1.0 

6.9 

0.86 

Servo 

10.0 

1.0 

6.9 

0.86 


Note 1: Baseline RQAS Frequency and Damping are based on 

roots of the W' -Plane. 


Optimal Gain Matrix : 


K 


1.8623 -1.2121 

0.4506 1.8538 


-0.7767 

-2.9208 


-1.3093 

1.4859 


Weighting Matrices — Diagonal Elements : 


Q' = (0.05, 10.0, 0.1* 0.75, 10.0) 
R' = (5.0, 1.0) 
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Table B.37 Summary of Lateral-Directional RQAS Response 
for a Sea-Level Climb — Output-Weighting 


Time Domain Response ; 

Open Baseline Percent 

Loo P RQAS Reduction 



Max. 

RMS 

Max. 

RMS 

Max. 

RMS 

a y 

(f/s2) 

6.1 

2.7 

3.4 

1.1 

44% 

59% 

B 

(deg) 

6.2 

2.8 

6.0 

1.9 

3% 

32% 

P 

(deg/s) 

6.4 

2.7 

1.8 

0.6 

72% 

78% 

r 

(deg/s) 

10.6 

4.8 

6.8 

2.7 

36% 

44% 

4> 

(deg) 

3.2 

1.3 

1.7 

0.7 

47% 

46% 

«df 

(deg) 

- 

- 

6.2 

2.8 

- 

- 

«sr 

(deg) 

— 

- 

14.7 

5.7 

- 

- 

, $df 
(deg/s) 

- 

- 

51.4 

- 

- 

- 

^sr 

(deg/s) 

- 

- 

27.1 

- 

- 

- 


B . 55 


Table B. 37 Summary of Lateral-Directional RQAS Response 
for a Sea-Level Climb — Output-Weighting 
(Continued) 


System Eigenvalues : 

Open Loop 
Frequency Damping 


Baseline RQAS 
Frequency Damping 


Spiral Mode 

0.02 

-1.0 

0.78 

1.0 

Roll Mode 

2.7 

1.0 

7.4 

1.0 

Dutch Roll 

2.1 

0.12 

2.2 

0 . 58 

Dutch Roll 

2.1 

0.12 

2.2 

0 . 58 

Servo 

10.0 

1.0 

6.8 

0.87 

Servo 

10.0 

1.0 

6.8 

0.87 

Note 1: Baseline RQAS 

roots of the 

Frequency 
W' -Plane . 

and Damping are 

based on 


Optimal Gain Matrix ; 



f— 



— 1 


2.0004 

-0.8556 

-0.6140 

-1.0563 

K = 

-0.0932 

1.1370 

-2.4060 

1.1638 


Weighting Matrices — Diagonal Elements : 

Q' * (0.05, 10.0, 0.1, 0.75, 10.0) 

R' = (7.0, 1.8) 
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Tab ■ e B. 38 Summary of Lateral-Directional RQAS Response 
for a Climb at 5000 ft — Output-Weighting 


Time Domain Response : 

Open Baseline Percent 

Loop RQAS Reduction 



Max. 

RMS 

Max. 

RMS 

Max. 

RMS 

a y 

<f/s2) 

2.9 

1.2 

1.9 

0.7 

34% 

42% 

B 

(deg) 

3.3 

1.3 

3.2 

1.2 

3% 

8% 

P 

(deg/s) 

2.6 

1.1 

0.7 

0.3 

73% 

73% 

r 

(deg/s) 

3.9 

1.9 

3.3 

1.4 

15% 

26% 

<t> 

(deg) 

1.1 

0.5 

0.6 

0.3 

45% 

40% 

, 6 df 

(deg) 

- 

- 

3.1 

1.5 

- 

- 

$sr 

(deg) 

- 

- 

9.3 

3.9 

- 

- 

, $df 
(deg/s) 

- 

- 

28.8 

- 

- 

- 

$sr 

(deg/s) 

- 

- 

22.8 

- 

- 

- 
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Table B.38 Summary of Lateral-Directional RQAS Response 
for a Climb at 5000 ft— Output -Weight ing 
(Continued) 


System Eigenvalues : 


Open Loop Baseline RQAS 


Frequency Damping Frequency Damping 


Spiral Mode 

0.02 

-1.0 

0.70 

1.0 

Roll Mode 

2.5 

1.0 

6*1 

1.0 

Dutch Roll 

2.1 

0.11 

2.4 

0.73 

Dutch Roll 

2.1 

0.11 

2.4 

0.73 

Servo 

10.0 

1.0 

6.8 

0.89 

Servo 

10.0 

1.0 

6.8 

0.89 

Note 1: Baseline RQAS 

roots of the 

Frequency 
W' -Plane. 

and Damping are 

based on 


Optimal Gain Matrix : 


^ — 

r* 


— 


1.7209 

-0.7476 

-0.6779 -0.8545 

K = 

0.1127 

1.5494 

-3.1532 1.4838 


Weighting Matrices — Diagonal Elements : 


Q' = (0.05, 10.0, 0.1, 0.75, 10.0) 
R' = (10.0, 1.0) 
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Table B.39 Summary of Lateral-Directional RQAS Response 
for Cruise at 20000 ft — Output-Weighting 


Time Domain Response: 


Open 

Loop 


Baseline 

RQAS 


Percent 

Reduction 



Max . 

RMS 

Max. 

RMS 

Max . 

RMS 

a y 

(f/s2) 

2.4 

1.1 

1.5 

0.4 

38% 

64% 

8 

(deg) 

1.9 

0.8 

1.6 

0.5 

16% 

38% 

P 

(deg/s) 

2.3 

1.0 

1.0 

0.4 

57% 

60% 

r 

(deg/s) 

4.1 

1.8 

1.8 

0.8 

56% 

56% 

♦ 

(deg) 

1.0 

0.4 

0.6 

0.3 

40% 

25% 

6 df 

(deg) 

- 

- 

2.2 

0.9 

- 

- 

^sr 

(deg) 

- 

- 

5.2 

2.1 

- 

- 

, $df 
(deg/s) 

- 

- 

18.2 

- 

- 

- 

$sr 

(deg/s) 

- 

- 

22.8 

- 

- 

- 
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Table B.39 Summary of Lateral-Directional RQAS Response 
for Cruise at 20000 ft — Output-Weighting 
(Continued) 


System Eigenvalues : 

Open Loop 
Frequency Damping 


Baseline RQAS 
Frequency Damping 


Spiral Mode 

0.003 

-1.0 

1.1 

1.0 

Roll Mode 

2.4 

1.0 

2.7 

1.0 

Dutch Roll 

2.6 

0.08 

3.9 

0.70 

Dutch Roll 

2.6 

0.08 

3.9 

0.70 

Servo 

10.0 

1.0 

8.0 

0.81 

Servo 

10.0 

1.0 

8.0 

0.81 


Note 1: Baseline RQAS Frequency and Damping are based on 

roots of the W' -Plane. 


Opt 

K 


imal Gain Matrix : 

2.6890 -0.6745 

-1.4528 1.2040 


-0.8628 

-2.7719 


-0.6222 

1.6245 


Weighting Matrices — Diagonal Elements : 


Q' * (0.05, 10.0, 0.1, 0.75, 10.0) 
R’ = (10.0, 1.0) 
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Table B.40 Summary of Lateral-Directional RQAS Response 

for an Approach at Sea-Level — Output-Weighting 

Time Domain Response: 


Open Baseline Percent 

Loop RQAS Reduction 



Max. 

RMS 

Max. 

RMS 

Max . 

RMS 

a y 

(f/s2) 

5.8 

2.4 

2.6 

1.1 

50% 

54% 

8 

(deg ) 

10.2 

4.2 

5.9 

2.6 

42% 

38% 

P . 

(deg/s) 

7.4 

3.1 

1.6 

0.7 

78% 

77% 

r 

(deg/s) 

10.5 

5.8 

6.6 

3.4 

37% 

41% 

♦ 

(deg) 

3.8 

2.0 

0.9 

0.4 

76% 

56% 

$df 

(deg) 

- 

- 

4.7 

2.1 

— 

— 

«sr 

(deg) 

- 

- 

12.1 

6.2 

— 

— 

«df 

(deg/s) 

- 

- 

17.8 

— 

— 

— 

fisr 

(deg/s) 

- 

- 

24.5 

- 

- 

- 
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Table B.40 Summary of Lateral-Directional RQAS Response 

for an Approach at Sea-Level — Output-Weighting 
(Continued) 


System Eigenvalues ; 



Open 

Loop 

Baseline 

RQAS 


Frequency Damping 

Frequency 

Damp i 

Spiral Mode 

0.003 

-1.0 

0.65 

1.0 

Roll Mode 

2.1 

1.0 

2.6 

1.0 

Dutch Roll 

1.6 

0.11 

1.7 

0.40 

Dutch Roll 

1.6 

0.11 

1.7 

0.40 

Servo 

10.0 

1.0 

8.5 

1.00 

Servo 

10.0 

1.0 

9.0 

1.00 


Note 1: Baseline RQAS Frequency and Damping are based on 

roots of the W' -Plane. 


Optimal Gain Matrix : 


K 


0.8230 -0.6195 -0.2863 -0.9670 

0.2765 0.4704 -1.8866 0.4482 


Weighting Matrices — Diagonal Elements : 


Q' = (0.05, 10.0, 0.1, 0.75, 10.0) 
R' = (10.0, 3.0) 
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Table B.41 Longitudinal Control-Rate-Weighting 

Fixed Gain Stability for Sea-Level Take-Off 


Baseline Eigenvalues: 


Mode 

Frequency 

Damping 

Phugoid 

0.0940 

0.868 

Phugoid 

0.0940 

0.868 

Short Period 

6.28 

0.700 

Short Period 

6.28 

0.700 

Servo 

2.38 

1.0 

Servo 

11.6 

1.0 

Filter 

14.6 

1.0 

Filter 

36.6 

1.0 


Fixed Gain Eigenvalues: 


Foward c.g. Mid c.g. Aft c.g. 


Freq. 

Damp. 

Freq. 

Damp. 

Fre q* 

Damp. 

0.222 

0.309 

0.224 

0.311 

0.221 

0.318 

0.222 

0.309 

0.224 

0.311 

0.221 

0.318 

6.51 

0.743 

6.13 

0.723 

5.81 

0.702 

6.51 

0.743 

6.13 

0.723 

5.81 

0.702 

2.40 

1.0 

2.71 

1.0 

3.03 

1.0 

11.5 

1.0 

11.5 

1.0 

11.5 

1.0 

15.2 

1.0 

15.0 

1.0 

14.9 

1.0 

63.6 

1.0 

63.8 

1.0 

63.9 

1.0 
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Table B.42 Longitudinal Control-Rate-Weighting 

Fixed Gain Stability for Sea-Level Climb 


Baseline Eigenvalues: 


Mode 

Frequency 

Damp i ng 

Phugoid 

0.0874 

0.791 

Phugoid 

0.0874 

0.791 

Short Period 

6.02 

0.682 

Short Period 

6.02 

0.682 

Servo 

3.47 

1.0 

Servo 

11.8 

1.0 

Filter 

15.3 

1.0 

Filter 

63.8 

1.0 


Fixed Gain Eigenvalues: 


Foward c.g. Mid c.g. Aft c.g. 


Freq^ 

Damp. 

Freq. 

Damp. 

Freq. 

Damp . 

0.213 

0.368 

0.214 

0.370 

0.214 

0.372 

0.213 

0.368 

0.214 

0.370 

0.214 

0.372 

6.80 

0.699 

6.44 

0.671 

6.14 

0.643 

6.80 

0.699 

6.44 

0.671 

6.14 

0.643 

2.56 

1.0 

2.88 

1.0 

3.17 

1.0 

11.7 

1.0 

11.7 

1.0 

11.7 

1.0 

15.9 

1.0 

15.7 

1.0 

15.8 

1.0 

63.5 

1.0 

63.7 

1.0 

63.9 

1.0 
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Table B.43 Longitudinal Control-Rate-Weighting 

Fixed Gain Stability for Climb at 5000 ft 


Baseline Eigenvalues: 


Mode 

Frequency 

Damping 

Phugoid 

0.0837 

0.798 

Phugoid 

0.0837 

0.798 

Short Period 

6.13 

0.633 

Short Period 

6.13 

0.633 

Servo 

3.48 

1.0 

Servo 

11.7 

1.0 

Filter 

15.6 

1.0 

Filter 

64.1 

1.0 


Fixed Gain Eigenvalues: 


Foward c.g. Mid c.g. Aft c.g. 


Freq. 

Damp. 

Freq. 

Damp. 

Freq. 

Damp. 

0.197 

0.372 

0.195 

0.378 

0.198 

0.376 

0.197 

0.372 

0.195 

0.378 

0.198 

0.376 

6.59 

0.675 

6.23 

0.642 

5.95 

0.610 

6.59 

0.675 

6.23 

0.642 

5.95 

0.610 

2.80 

1.0 

3.14 

1.0 

3.46 

1.0 

11.6 

1.0 

11.6 

1.0 

11.6 

1.0 

15.6 

1.0 

15.5 

1.0 

15.4 

1.0 

63.5 

1.0 

63.7 

1.0 

63.9 

1.0 
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Table B . 44 Longitudinal Control-Rate-Weighting 

Fixed Gain Stability for Cruise at 20,000 ft 


Baseline Eigenvalues: 


Mode 

Frequency 

Damping 

Phugoid 

0.0806 

0.778 

Phugoid 

0.0806 

0.778 

Short Period 

10.2 

0.481 

Short Period 

10.2 

0.481 

Servo 

1.27 

1.0 

Servo 

11.5 

1.0 

Filter 

27.9 

1.0 

Filter 

126 

1.0 


Fixed Gain Eigenvalues: 


Foward c.g. Mid c.g. Aft c.g. 


Freq. 

Damp. 

Freq. 

Damp. 

Freq. 

Damp. 

0.143 

0.513 

0.143 

0.517 

0.146 

0.508 

0.143 

0.513 

0.143 

0.517 

0.146 

0.508 

6.97 

0.510 

6.69 

0.461 

6.48 

0.417 

6.97 

0.510 

6.69 

0.461 

6.48 

0.417 

3.48 

1.0 

3.79 

1.0 

4.03 

1.0 

11.6 

1.0 

11.6 

1.0 

11.6 

1.0 

16.5 

1.0 

16.4 

1.0 

16.3 

1.0 

63.4 

1.0 

63.7 

1.0 

63.9 

1.0 
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Table B.45 Longitudinal Control-Rate-Weighting 

Fixed Gain Stability for Sea-Level Approach 


Baseline Eigenvalues: 


Mode 

Frequency 

Damping 

Phugoid 

0.187 

0.934 

Phugoid 

0.187 

0.934 

Short Period 

7.36 

0.667 

Short Period 

7.36 

0.667 

Servo 

3.22 

1.0 

Servo 

10.6 

1.0 

Filter 

21.4 

1.0 

Filter 

32.4 

1.0 


Fixed Gain Eigenvalues: 


Foward 

c.g. 

Mid 

c.g. 

Aft 

c.g. 

Freq. 

Damp. 

Freq. 

Damp. 

Freq. 

Damp. 

0.286 

0.518 

0.289 

0.582 

0.294 

0.648 

0.286 

0.518 

0.289 

0.582 

0.294 

0.648 

5.56 

0.727 

5.24 

0.745 

4.97 

0.767 

5.56 

0.727 

5.24 

0.745 

4.97 

0.767 

3.97 

1.0 

3.86 

1.0 

3.71 

1.0 

10.4 

1.0 

10.4 

1.0 

10.4 

1.0 

14.2 

1.0 

14.1 

1.0 

14.0 

1.0 

64.1 

1.0 

64.1 

1.0 

64.2 

1.0 
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Table B. 46 Lateral Directional Control-Rate-Weighting 
Fixed Gain Stability for Sea-Level Take-Off 



Baseline 

Fixed 

Gain 


Frequency 

Damping 

Frequency 

Damping 

Spiral 

0.361 

1.0 

0.485 

1.0 

Roll 

2.76 

1.0 

4.14 

1.0 

Dutch Roll 

2.04 

0.550 

2.07 

0.558 

Dutch Roll 

2.04 

0.550 

2.07 

0.558 

Servo 

5.87 

0.855 

5.01 

0.972 

Servo 

5.87 

0.855 

5.01 

0.972 

Filter 

11.2 

1.0 

12.1 

1.0 

Filter 

25.5 

1.0 

21.4 

1.0 


Table B.47 Lateral Directional Control-Rate-Weighting 
Fixed Gain Stability for Sea-Level Climb 



Baseline 

Fixed 

Gain 


Frequency 

Damping 

Frequency 

Damping 

Spiral 

0.594 

1.0 

0.570 

1.0 

Roll 

3.25 

1.0 

3.06 

1.0 

Dutch Roll 

2.75 

0.514 

2.77 

0.507 

Dutch Roll 

2.75 

0.514 

2.77 

0.507 

Servo 

5.47 

0.882 

5.59 

0.871 

Servo 

5.47 

0.882 

5.59 

0.871 

Filter 

12.6 

1.0 

12.4 

1.0 

Filter 

21.3 

1.0 

21.8 

1.0 
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Table B.48 Lateral Directional Control-Rate-Weighting 
Fixed Gain Stability for Climb at 5000 ft 



Baseline 

Fixed 

Gain 


Frequency 

Damping 

Frequency 

Damp i ng 

Spiral 

0.581 

1.0 

0.613 

1.0 

Roll 

3.19 

1.0 

3.11 

1.0 

Dutch Roll 

2.72 

0.503 

2.74 

0.498 

Dutch Roll 

2.72 

0.503 

2.74 

0.498 

Servo 

5.57 

0.847 

5.41 

0.880 

Servo 

5.57 

0.847 

5.41 

0.880 

Filter 

12.6 

1.0 

12.4 

1.0 

Filter 

21.8 

1.0 

21.8 

1.0 




Table B . 49 

i 

! 


Lateral Directional Control-Rate-Weighting 
Fixed Gain Stability for Cruise at 20,000 ft 


I 



Baseline 

Fixed 

Gain 


Frequency 

Damp i ng 

Frequency 

Damp i ng 

Spiral 

0.799 

1.0 

0.814 

1.0 

Roll 

2.55 

1.0 

2.39 

1.0 

Dutch Roll 

3.48 

0.449 

3.92 

0.333 

Dutch Roll 

3.48 

0.449 

3.92 

0.333 

Servo 

6.27 

0.716 

5.93 

0.717 

Servo 

6.27 

0.716 

5.93 

0.717 

Filter 

14.1 

1.0 

13.2 

1.0 

Filter 

23.6 

1.0 

22.6 

1.0 
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Table B.50 Lateral Directional Control-Rate-Weighting 
Fixed Gain Stability for Sea-Level Approach 



Baseline 

Fixed 

Gain 


Frequency 

Damping 

Frequency 

Damping 

Spiral 

0.546 

1.0 

0.511 

1.0 

Roll 

2.74 

1.0 

3.31 

1.0 

Dutch Roll 

1.85 

0.404 

1.75 

0.468 

Dutch Roll 

1.85 

0.404 

1.75 

0.468 

Servo 

5.12 

1.0 

4.40 

1.0 

Servo 

7.61 

1.0 

7.03 

1.0 

Filter 

11.9 

1.0 

11.7 

1.0 

Filter 

19.4 

1.0 

21.1 

1.0 
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Table B.51 Longitudinal Output-Weighting Fixed 

Gain Stability for Sea-Level Take-Off 


Baseline Eigenvalues: 


Mode 

Frequency 

Damping 

Phugoid 

0.0895 

1.0 

Phugoid 

0.111 

1.0 

Short Period 

8.72 

0.866 

Short Period 

8.72 

0.866 

Servo 

2.09 

1.0 

Servo 

11.1 

1.0 


Fixed Gain Eigenvalues: 


Foward 

c.g. 

Mid 

c.g. 

Aft 

c.g. 

Freq. 

Damp. 

Freq. 

Damp. 

Freq. 

Damp. 

0.224 

0.350 

0.226 

0.350 

0.224 

0.355 

0.224 

0.350 

0.226 

0.350 

0.224 

0.355 

11.2 

0.704 

10.8 

0.695 

10.5 

0.689 

11.2 

0.704 

10.8 

0.695 

10.5 

0.689 

1.70 

1.0 

1.82 

1.0 

1.92 

1.0 

11.1 

1.0 

11.1 

1.0 

11.2 

1.0 




B . 71 


Table B.52 Longitudinal Output-Weighting Fixed 
Gain Stability for Sea-Level Climb 


Baseline Eigenvalues: 


Mode 

Frequency 

Damping 

Phugoid 

0.0887 

0.886 

Phugoid 

0.0887 

0.886 

Short Period 

11.6 

0.661 

Short Period 

11.6 

0.661 

Servo 

1.94 

1.0 

Servo 

11.3 

1.0 


Fixed Gain Eigenvalues: 


Foward c.g. Mid c.g. Aft c.g. 


Freq. 

Damp. 

Freq. 

Damp. 

Freq. 

Damp. 

0.216 

0.420 

0.216 

0.419 

0.217 

0.418 

0.216 

0.420 

0.216 

0.419 

0.217 

0.418 

12.1 

0.668 

11.7 

0.658 

11.4 

0.649 

12.1 

0.668 

11.7 

0.658 

11.4 

0.649 

1.77 

1.0 

1.89 

1.0 

1.99 

1.0 

11.3 

1.0 

11.3 

1.0 

11.3 

1.0 
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Table B.53 Longitudinal Output-Weighting Fixed 
Gain Stability for Climb at 5000 ft 


Baseline Eigenvalues: 


Mode 

Frequency 

Damping 

Phugoid 

0.0849 

0.885 

Phugoid 

0.0849 

0.885 

Short Period 

11.5 

0.636 

Short Period 

11.5 

0.636 

Servo 

2.02 

1.0 

Servo 

11.2 

1.0 


Fixed Gain Eigenvalues: 


Foward 

c.g. 

Mid 

c.g. 

Aft 

c.g. 

Freq. 

Damp . 

Freq. 

Damp. 

Freq. 

Damp. 

0.199 

0.424 

0.198 

0.428 

0.201 

0.422 

0.199 

0.424 

0.198 

0.428 

0.201 

0.422 

11.8 

0.653 

11.5 

0.643 

11.2 

0.634 

11.8 

0.653 

11.5 

0.643 

11.2 

0.634 

1.86 

1.0 

1.98 

1.0 

2.09 

1.0 

11.2 

1.0 

11.2 

1.0 

11.3 

1.0 
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Table B.54 Longitudinal Output-Weighting Fixed 

Gain Stability for Cruise at 20,000 ft 


Baseline Eigenvalues: 


Mode 

Frequency 

Damping 

Phugoid 

0.0772 

0.809 

Phugoid 

0.0772 

0.809 

Short Period 

13.6 

0.514 

Short Period 

13.6 

0.514 

Servo 

1.93 

1.0 

Servo 

11.3 

1.0 


Fixed Gain Eigenvalues: 


Foward 

c.g. 

Mid 

c.g. 

Aft 

c.g. 

Freq. 

Damp. 

Freq. 

Damp. 

Freq. 

Damp. 

0.145 

0.585 

0.145 

0.585 

0.148 

0.571 

0.145 

0.585 

0.145 

0.585 

0.148 

0.571 

13.3 

0.542 

12.9 

0.529 

12.6 

0.518 

13.3 

0.542 

12.9 

0.529 

12.6 

0.518 

2.17 

1.0 

2.30 

1.0 

2.41 

1.0 

11.2 

1.0 

11.2 

1.0 

11.2 

1.0 
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Table B.55 Longitudinal Output-Weighting Fixed 
j Gain Stability for Sea-Level Approach 

i 

I 


Baseline Eigenvalues: 


Mode 

Frequency 

Damp i ng 

Phugoid 

0.178 

0.939 

Phugoid 

0.178 

0.939 

Short Period 

6.89 

0.861 

Short Period 

6.89 

0.861 

Servo 

4.50 

1.0 

Servo 

10.3 

1.0 


Fixed Gain Eigenvalues: 


Foward c.g. Mid c.g. Aft c.g. 


Freq. 

Damp. 

Freq. 

Damp. 

Freq. 

Damp. 

0.281 

0.543 

0.283 

0.598 

0.287 

0.651 


0.281 

0.543 

0.283 

0.598 

0.287 

0.651 

9.57 

0.731 

9.27 

0.739 

9.03 

0.746 

9.57 

0.731 

9.27 

0.739 

9.03 

0.746 

2.61 

1.0 

2.45 

1.0 

2.30 

1.0 

10.3 

1.0 

10.3 

1.0 

10.3 

1.0 
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Table B.56 Lateral Directional Output-Weighting 

Fixed Gain Stability for Sea-Level Take-Off 



Baseline 

Fixed 

Gain 


Frequency 

Damping 

Frequency 

Damping 

Spiral 

0.669 

1.0 

0.643 

1.0 

Roll 

7.47 

1.0 

7.89 

1.0 

Dutch Roll 

1.95 

0.606 

1.77 

0.561 

Dutch Roll 

1.95 

0.606 

1.77 

0.561 

Servo 

6.90 

0.862 

6.58 

0.908 

Servo 

6.90 

0.862 

6.58 

0.908 


Table B.57 Lateral Directional Output-Weighting 

Fixed Gain Stability for Sea-Level Climb 



Baseline 

Fixed 

Gain 


Frequency 

Damping 

Frequency 

Damping 

Spiral 

0.776 

1.0 

0.683 

1.0 

Roll 

7.39 

1.0 

7.06 

1.0 

Dutch Roll 

2.22 

0.578 

2.30 

0.586 

Dutch Roll 

2.22 

0.578 

2.30 

0.586 

Servo 

6.84 

0.874 

7.31 

0.833 

Servo 

6.84 

0.874 

7.31 

0.833 
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Table B . 58 Lateral Directional Output-Weighting 

Fixed Gain Stability for Climb at 5000 ft 



Baseline 

Fixed 

Gain 


Frequency 

Damping 

Frequency 

Damping 

Spiral 

0.699 

1.0 

0.715 

1.0 

Roll 

6.08 

1.0 

7.06 

1.0 

Dutch Roll 

2.44 

0.728 

2.30 

0.572 

Dutch Roll 

2.44 

0.728 

2.30 

0.572 

Servo 

6.78 

0.890 

7.17 

0.837 

Servo 

6.78 

0.890 

7.17 

0.837 


Table B. 59 Lateral Directional Output-Weighting 

Fixed Gain Stability for Cruise at 20,000 ft 



Baseline 

Fixed 

Gain 


Frequency 

Damping 

Frequency 

Damping 

Spiral 

1.09 

1.0 

0.822 

1.0 

Roll 

2.67 

1.0 

4.61 

1.0 

Dutch Roll 

3.92 

0.702 

3.82 

0.606 

Dutch Roll 

3.92 

0.702 

3.82 

0.606 

Servo 

7.97 

0.808 

8.28 

0.722 

Servo 

7.97 

0.808 

8.28 

0.722 
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Table B.60 Lateral Directional Output-Weighting 

Fixed Gain Stability for Sea-Level Approach 



Baseline 

Fixed 

Gain 


Frequency 

Damping 

Frequency 

Damping 

Spiral 

0.652 

1.0 

0.643 

1.0 

Roll 

2.63 

1.0 

5.58 

1.0 

Dutch Roll 

1.68 

0.403 

1.62 

0.467 

Dutch Roll 

1.68 

0.403 

1.62 

0.467 

Servo 

8.49 

1.0 

6.01 

1.0 

Servo 

9.06 

1.0 

8.34 

1.0 
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APPENDIX C. 


MODIFICATION DRAWINGS 


This appendix contains a full set of preliminary 
structural drawings for the proposed Cessna 402B modifi- 
cations discussed in Chapter 4. Note that all of these 
drawings have been reduced from the originals and there- 
fore, the indicated scales no longer apply. The origi- 
nal drawings are on file at the University of Kansas 
Flight Research Laboratory. 
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